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Abstract 

Constructing a well-posed variational principle is a non-trivial issue in general relativity. For 
spacelike and timelike boundaries, one knows that the addition of the Gibbons-Hawking-York (GHY) 
counter-term will make the variational principle well-defined. This result, however, does not directly 
generalize to null boundaries on which the 3-metric becomes degenerate. In this work, we address 
the following question: What is the counter-term that may be added on a null boundary to make 
the variational principle well-defined? We propose the boundary integral of 2^—g (0 -|- k) as an ap¬ 
propriate counter-term for a null boundary. We also conduct a preliminary analysis of the variations 
of the metric on the null boundary and conclude that isolating the degrees of freedom that may be 
fixed for a well-posed variational principle requires a deeper investigation. 


1 Introduction and Summary 

Just like any other field theory, the dynamics of gravity can be obtained from an action, the Einstein- 
Hilbert action. On varying the action with respect to the metric, we obtain an equations-of-motion term 
and a boundary term generated by integration by parts. The equations of motion turn out to be second 
order in the derivatives of the metric. But the boundary term is unusual, as it contains variations both 
of the metric and its first derivatives, both tangential and normal to the boundary [1]. Thus, setting the 
variation of the action to zero will lead to the equations of motion on the bulk only if we fix the metric 
and its normal derivatives on the boundary. (The tangential derivatives get fixed when the metric is 
fixed on the boundary.) The problem with such a structure is that it makes the variational principle 
ill-defined [2]. In general, the equations of motion and the boundary conditions will turn out to be 
inconsistent [2, 3]. 

There is a widely accepted prescription for resolving this issue. One adds an extra boundary term to 
the action, called a counter-term, such that the surface term in the variation of the new action contains 
only variations of the metric and does not involve the variations of the normal derivatives. (We shall be 
using the terms surface term, boundary term and counter-term interchangeably. Note that the use of 
the word “counter-term” here, unlike in many places in the literature, does not mean that our boundary 
term is supposed to cancel off a divergent part of the action.) Thus, we need to fix only the metric on the 
boundary and the variational principle becomes well-posed. The most commonly used counter-term is 
the Gibbons-Hawking-York (GHY) counter-term [4, 5], although there are other counter-terms available 
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in the literature [6]. In cases where the boundary is taken to infinity, other boundary terms in addition 
to the Gibbons-Hawking-York term are added to make the action finite on classical solutions [7, 8]. But 
we shall not be discussing this issue in the current work. 

The GHY prescription has been around for a long time and is now textbook material (see, e.g., [1, 7]). 
If the spacetime region V under consideration has a boundary surface denoted by 9V, the counter-term 
is an integral over dV of essentially the product of the square root of the determinant of the 3-metric 
on the surface and the divergence of the unit normal to the surface. The procedure of its construction 
uses (a) the unit normal to dV, (b) the induced metric hab on dV and (c) the covariant derivative on 
dV compatible with the induced metric hat- These structures are well-defined for spacelike or timelike 
regions of the boundary surface. The trouble is, we can have a third kind of region on the boundary 
surface, viz. a null surface, which is ubiquitous in general relativity. The horizons of black holes, for 
example, are null surfaces. It is generally accepted that entropy and temperature can be associated 
with black hole horizons [9, 10, 11, 12, 13]. There is also a claim in the literature that entropy and 
temperature can also be ascribed to any local Rindler horizon [14, 15], which is also a null surface. These 
properties, as far as we know, are not shared by any spacelike or timelike surface. Some physical cases 
where we might require a counter-term for a null surface are the case of a spacetime with a black hole 
where we want to consider the black hole horizon as one of the boundaries, the case of the interior of 
a causal diamond formed by the intersection of the future light cone from a point P and the past light 
cone from a point Q in the future of P [16], etc. The trouble with the standard prescription is that it is 
not directly applicable to a null surface as the normal to a null surface has a zero norm and the natural 
3-metric on a null surface is degenerate. Also, the notion of a covariant derivative on a non-null surface 
does not naturally extend to a null surface, making the usual procedure inapplicable. 

A possible work-around is to consider a timelike surface infinitesimally separated from the null surface 
(“stretched horizon”), perform the calculations and then take the null limit. As we shall demonstrate, this 
approach does give a proper, finite expression for the boundary term and the counter-term when the null 
limit is taken with some reasonable assumptions about the limiting behaviour of the metric components. 
But, given the importance of this problem, one would like to have a first-principle approach, leading to 
a prescription for the counter-term to be added on the null surface purely based only on properties of 
the surface, without recourse to any limiting procedure. Surprisingly, this issue does not seems to have 
attracted sufficient attention in the literature. As far as we know, there is no approach available in the 
literature that allows one to start from first principles and find the counter-term for a null surface. (The 
only related reference we could find was [17] which treats a very specialized metric and is not applicable 
to a general null surface and does not address several important technical and conceptual issues.) 

There is another important reason for attempting this study. For the non-null surfaces, we also obtain, 
from the variational principle itself, the result that only the components of the induced metric (modulo 
diffeomorphisms) need to be fixed at the boundary [18]. So, in a timelike/spacelike surface one can take 
the gravitational degrees of freedom as being contained in the 6 components of the induced 3-metric hab- 
Given the freedom in transforming the 3 coordinates on dV, this leaves 3 independent degrees of freedom. 
(It is possible to reduce it further to 2 using the freedom in foliation, giving two degrees of freedom for 
the gravitational field per point [19, 20].) When we try to obtain the results for a null surface as a limit 
of the results on a non-null surface, the question arises as to how to characterize these degrees of freedom 
in the null limit. The properly defined induced metric for null surfaces is essentially a 2-metric qab- It is 
not clear how the six degrees of freedom in the induced metric hab on the non-null surface could relate 
to the three degrees of freedom in the 2-metric Qab on taking the null limit. 

We stress that it should be possible to tackle the above issues from first principles, by just studying 
the form of the variation of the action 5Ag^ on the null boundary. (We shall use the notation SAgy for 
the boundary term in the variation of the Einstein-Hilbert action for a boundary dV of the spacetime 
region V under consideration.). The algebraic structure of SAgy should suggest (a) the counter-term 
that may be added to make the variation well-defined and (b) the variables which need to be fixed on the 
boundary. (For the corresponding analysis in the case of a non-null boundary, see [18].) It is important 
to perform this computation for a null boundary. 

We should also point out that a complete analysis of the well-posedness of the action requires con¬ 
sideration of the equations of motion as well, since well-posedness essentially requires that the boundary 
conditions be consistent with the equations of motion and that each choice of boundary conditions selects 
out a unique, different solution of the equations of motion [2]. The evolution of geometry from one spatial 
surface to another under the Einstein equations is well-understood and there is a standard, accepted for¬ 
malism in the Arnowitt-Deser-Misner (ADM) formalism [21]. On the other hand, the evolution of geome¬ 
try starting from a null surface has been analyzed by several authors [22, 23, 24, 25, 26, 27, 28, 29, 30, 31], 
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but no one approach has yet come to be accepted as the standard. Details like the specification of the 
degrees of freedom and the identification of constraints vary from approach to approach. One of these 
formalisms has to be chosen and then the well-posedness of the action has to be analyzed in its framework. 
We intend to pursue this matter in a future publication. 

We will provide comprehensive discussion of these and related issues from several different perspectives 
in this paper. Because of the rather extensive nature of the discussion, we hrst summarize the key results 
obtained in this work. 


• We briefly review the procedure of obtaining a counter-term and identifying the variables which 
need to be fixed on the boundary just from examining the algebraic structure of SAgy on non-null 
boundary. Then, we perform a very similar analysis and obtain 5Agy on a general null surface. 
The final result is given by 

^AnuU = [ X {—26 [^/—g{& + k)] + 'f—g [©ab — (0 -f K)qab\ 6q°'^ + Pc6i'^{ . (1) 

JdV 

(There is also an ignorable total derivative term on the 3-surface, which has been omitted here for 
simplicity but can be found in the main discussion later on.) The general result as given in Eq. (1) 
can be interpreted by the following identifications: ^a = daf' is the normal to the </> = constant 
boundary surface which is null (other surfaces of constant (j) may not be null), is an auxiliary null 
vector and qab = gab + (-a^b + kaib is the induced metric. The symmetric tensor &ab = q^qb^m^n 
is the second fundamental form, 0 = 0“ is the expansion scalar and k is the non-affinity coefficient 
on the null surface, i.e. (Using ^a = da4>, as opposed to ia = Adaf) for some scalar 

A, will render k = 0 if all the (j = constant surfaces are null but not when (j) = constant is null 
only for a specific value of (f; this is why k 0 even though ia = dafi; see discussion later.). In 
a coordinate system adapted to a null surface with a particular choice of ka, we can write this 
decomposition in the form 


[SAv]null= [ dXd^Z^{-2d[^{Q + K)] + ^[Qab-{Q + K)qab]6q‘^'’ + Pc6P} , (2) 

JdV 

The integration is over a parameter A, varying along the null geodesic congruence such that ka = 
VaA locally, and two coordinates that are constant along the null geodesics, with d‘^z± = 

dz^dz'^. The symbol q is used for the determinant of the 2-metric qAB- The conjugate momentum 
to is y/q [0ab — (0 -f K)qab] and the conjugate momenta to turns out to be 

Pc = -VQkb[^ci’’ + ^^ic-2S^aiQ + f^)] ■ (3) 

We have also carried out the analysis for ia having the form Adaf) in Appendix G. In that case 
also, we obtain a decomposition of the same structure as above. 

The structure of [dAgyj^y^n is similar to the familiar decomposition that we obtain for the surface 
term on a non-null surface: 


[5AI, 


dV inon—null 



S{2^\K) - ^MiKab - Khab) Sh‘^^ 


(4) 


(Again, there is a total 3-derivative on the surface that we have omitted here but can be found 
in the main text.) Thus, Eq. (I) suggests that the counter-term to be added on a null surface is 
2y^(0 -I- k) {2y^[Q + k] for Eq. (2)) and that the intrinsic quantities to be fixed on the null 
surface are q°‘^ and i^". 

• We verify the result in Eq. (2) by deriving it through a different route. We start from the standard 
decomposition of the surface term for a non-null boundary and then take a careful null limit with 
some reasonable limiting behaviour for the metric components. We obtain back Eq. (2), reaffirming 
our faith in the result. 


• In an appendix, we illustrate these general results in two natural coordinate systems adapted to 
describe the metric near an arbitrary null surface. We first review the explicit construction of both 
these coordinate systems, which we call the Gaussian Null Coordinates (GNC) and Null Surface 
Foliation (NSF). The GNC metric has the feature that the null surface arises naturally as the limit 
of a sequence of non-null surfaces. Thus, we can easily apply the GHY prescription on a nearby 
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non-null surface and take the null limit. In the NSF, the fiducial null surface appears as one member 
of a set of null surfaces. Hence, the limiting procedure cannot be applied directly (and would need 
the explicit construction of an infinitesimally separated non-null surface). But the NSF has some 
other advantages. In both cases, we find that we obtain explicit results in agreement with Eq. (2). 

• The question of whether the variational principle is well-posed with this counter-term requires 
deeper analysis. Also, the limitations of our framework, since we impose certain constraints on the 
variations of the metric to keep, for example, the null surface as a null surface, does not allow us 
to directly connect these boundary conditions with the degrees of freedom in the theory. 

The conventions used in this paper are as follows: We use the metric signature (—, -I-,-b,-I-). The 
fundamental constants G, h and c have been set to unity. The Latin indices, a, 6,..., run over all space- 
time indices, and are hence summed over four values. Greek indices, a, f3 ,are used when we specialize 
to indices corresponding to a codimension-1 surface, i.e a 3—surface, and are summed over three values. 
Upper case Latin symbols. A, B ,..are used for indices corresponding to two-dimensional hypersurfaces, 
leading to sums going over two values. We shall use := to indicate equalities that are valid only on the 
null surface. 


2 The information contained in the boundary term of the vari¬ 
ation of the action 

The variation of an action, generically, contains two terms: (a) A bulk term, the vanishing of which will 
give us the equations of motion and (b) a boundary term, which is usually ignored. But the boundary 
term contains significant amount of information! By a careful analysis of this term, we can determine 
the structure of the counter-term that may be added to the action principle to make it well-defined and, 
with that counter-term, what degrees of freedom are to be fixed at the boundary. Some of these degrees 
of freedom may be eliminated by using gauge freedom and, once the gauge freedom is exhausted, the 
number of the remaining degrees of freedom represents the number of true degrees of freedom in the 
theory. Let us illustrate this idea in electromagnetism and gravity. 

2.1 Warm-up: Electromagnetic field 

As a warm up, let us consider the simple case of electromagnetism (see [19], Chapter 21). The action 
for the free electromagnetic field is given by 

A = - ^ / F^kF^^d^x- F,k = d,Ak - dkA, (5) 

lOTT Jy 

Let us begin by assuming that Ai are the dynamical variables and vary them in the action. This leads 
to 


SA=-^ [ IdkF^'^] SA,d^x + ^ [ F^HAdcTk = -^ [ \dkF^^]SA,d^x+^ / E • (6) 

^TtJv ItT Jgy ItT 7y L ^ AtT J 

where we have assumed that the boundary contributions arise from t = constant surfaces. Since bulk 
and boundary variations are independent, 5A = 0 will require vanishing of these two terms individually. 
Let us first assume (rather naively) that 5Ai is completely arbitrary in the bulk V. Then the vanishing 
of the bulk term will lead to the equations of motion dkF'^^ = 0. Once this is granted, we next want the 
on-shell boundary term to vanish. This is, of course, possible if we specify Ai on the boundary; but that 
is not required. The algebraic structure of the boundary term already tells us that we need not hx Aq at 
all. Further, even as regards A we only need to fix it modulo the addition of a gradient: A —^ A -b VF. 
This is because, E • VF = V • (EF) on-shell (since V • E = 0) which allows this term to be converted 
into a surface integral at spatial infinity and ignored. Thus, we only need to fix the part of A which is 
unaffected by the addition of a gradient; viz., the B = V x A. The magnetic field, since it is a spatial 
vector, has three components, of which one component is eliminated by the constraint V.B = 0 arising 
from the equations of motion. Thus, there are two degrees of freedom to be fixed at a boundary. Further, 
it can be verified that each choice of these two degrees of freedom selects out a single solution of the 
equations of motion. Thus, the action is well-posed and does not require the addition of a counter-term. 
In other words: 


4 


We do not tell the action principle what needs to he fixed at the boundary; instead, the action principle 
tells us what should be fixed. If the action is well-posed, the number of degrees of freedom that needs to 
he fixed at a boundary will correspond to the number of true degrees of freedom in the theory. 

In the case of the electromagnetism, it is the magnetic field that is to be fixed at the boundary. We 
naively thought it is Ai but the algebraic structure of the boundary term in the variation of the action 
corrects our mistake and tells us that it is the magnetic field. 


2.2 Gravity: Non-null boundary 

One can do exactly the same in the case of gravity. We will first assume (again naively) that gab are the 
dynamical variables and vary them in the Einstein-Hilbert action 


IbTrylEH = f R, (7) 

Jv 

where V indicates the spacetime region under consideration, over which the integration is performed. On 
varying the metric, we get 


IOttMeh = - [ d‘^x^/^G°‘''5gab + 

[ d^xy/^g'^^SR, 

Jv 

Jv 

= - [ d'^Xyd^G'^^Sgab + 

[ d'^x^/^ Vcw'^ 

Jv 

Jv 

= - [ d^x^/^ G°-^Sgab + 

f dfSucw'^ 

Jv 

JdV 

= - [ d'^XyG^ G°-^Sgab + 

f d^SQluc] 

Jv 

JdV 


where df S is the appropriate integration measure on the boundary, Uc is the normal with constant norm 
(0 or ±1) on the surface, and we have defined the function Q[Af\ for any 

one-form Ac as 


Q[Ac]=Ac{g'^HTl,-g^HTl^) . (9) 

which will prove to be useful throughout the paper. Straightforward algebra now leads to the following 
alternate expression for Q[nc]. 

Q[nc] = ^a{Snl) - (5(2V„n“) + VaUb 5g<^^- dnl = Sn'^ + g<^^8nb (10) 

This expression can be easily verified by expanding out and simplifying the RHS (see Appendix B for 
a derivation). Note that Snfua = 0, which implies that the vector 6nf lives on the boundary. The 
description so far is completely general and holds for both null and non-null boundaries with suitable 
measure UcdfS for integrating vector fields. 

For the non-null case ([18], also see Appendix B), we have the relation VaV°‘ = DaV^ — eabV^ valid 
for any which satisfies V^Ua = 0, where = n^'S/jUi. Further, the integration measure is now 
d^S = y/Tid^x with ria normalized to ±1 so that the surface term becomes the integral of d?x y/\h\Q[nc\. 
This will lead to the expression: 

Q[nc] = T>a(<5n“ ) - (5(2Von“) -t- (Van?, - euaab) Sh‘^’’ (11) 

Note that this algebraic manipulation has naturally led us to the combination {SIa^b — enaUb) without 
us having to introduce any geometrical considerations. This combination, of course, is the negative 
of the extrinsic curvature [1] defined by Kab = —hfVcnb and has the following properties: {a)Kij = 
Kjii {b)Kijn^ = 0; {c)K = g^^Kij = —Va'rA . The second property tells us that KijSg^^ = Kijdd^ which 
has been used to arrive at Eq. (11). The surface term is an integral of Q[nc] over df‘x^/\h\ which allows 
the first term Da{6nf) in Eq. (11) to be converted to a 2-dimensional surface integral and ignored. The 
rest of the terms will reduce to the following expression for the boundary term of EH action: 




non—null 



6{2^/\h\K) - VWliKab - Khab) 


( 12 ) 
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We need to be set to zero with appropriate boundary conditions to obtain the equations 

of motion in the bulk from the action principle 5A = 0. Unlike in the case of e.g., electromagnetism, 
this is not straightforward because we will need to set both 5[K\/h\ = 0 and 5h°‘^ = 0 on the boundary. 
This is not acceptable, since 5K contains variations of normal derivatives of the metric as well. Setting 
both 5[K\/h] = 0 and 5h°'^ = 0, in general, will lead to inconsistency with the Einstein’s equations. (For 
a discussion of this and other issues with fixing the metric and its derivatives at the boundaries, see, e.g, 
[ 2 ].) 

Incredibly enough, the algebraic structure of Eg. (12) itself suggests a way out. Since the term in¬ 
volving 5\K^/h] is the variation of a boundary integral, we can add this as a counter-term to the action 
such that the surface term arising from the new action does not depend on the variations of the normal 
derivatives of the metric. Once this is done, we only have to set Sh°‘^ = 0 on the boundary and this leads 
to a well-posed action principle. This leads to the GHY counter-term [4, 5], viz. integral of —2K\/h over 
the boundary, which is the most commonly used one for this purpose. 

Thus, the algebraic structure of the boundary term in the variation of the action itself can tell us not 
only (i) what counter-term can be added to the action to make the variation well-defined but also (ii) 
what degrees of freedom need to be fixed on the boundary [viz. hab] after the addition of this counter¬ 
term. Even here, eliminating four degrees of freedom using the gauge freedom that is the freedom of 
diffeomorphisms, one should end up with two remaining degrees [19], matching the number of degrees 
of freedom of the gravitational field obtained from considering gravitational radiation in the linearized 
approximation. 


3 Boundary Term of the Einstein-Hilbert Action: Null Surfaces 

The crucial difference between non-null and null boundaries is the following: If the surface is non-null, we 
can define an induced metric hab = gab — eriaUb and a covariant derivative operator Da compatible with 
it (i.e., Dahbc = 0). We can then re-express Va(<5n“ ) in terms of the covariant derivative Da{5n\) on 
the boundary which leads to significant simplification. We cannot do this for a null surface, which makes 
a vital algebraic difference. We, therefore, need to use some other procedure to analyse the algebraic 
structure of the boundary variation on a null-surface. 

One possible approach to handle a null surface would be to treat the null surface as the limit of 
a sequence of non-null surfaces, apply the prescription given above on the non-null surfaces and then 
take the limit in which this sequence goes over to the null surface. This can be done with the help of 
a parameter which labels the surfaces such that a particular value of the parameter corresponds to the 
null surface; we will do this in Section 3.6. However, we would also like to formulate a first-principle 
prescription that can be applied to the null surface as a surface in its own right and not as a fringe 
member of some family. This will be done in Section 3.3. As a preamble to these two derivations, we 
shall first recall several features of null surfaces relevant to our task. 

3.1 Null surfaces 

A null surface in a four-dimensional spacetime M is a three-dimensional submanifold N with the criterion 
that the metric obtained by the restriction of the full metric gab to the hypersurface N is degenerate, 
which means that it is possible to find non-zero vectors on N such that = 0. We shall assume 

that the full gab is non-degenerate, i.e every non-zero four-dimensional vector is mapped by the metric 
to a non-zero one-form. If v°‘ is the four-dimensional push-forward of and w'^ is the four-dimensional 
push-forward of w^, another vector at the same point on TV, then gabV°‘w^ = = 0. Thus, 

is orthogonal to every vector at that point on the surface including itself. A set of such null vectors on 
the surface, one vector per point, gives a null vector field normal to N. (There will be only one null 
curve on the surface passing through each point on A^. So the only freedom we have in the choice of the 
null vector is the freedom of scaling by a scalar factor.) Note that £“ has the peculiar property that it 
is both tangent and normal to the surface: tangent since we had started with the vector being on the 
null surface and hence we will be moving along the surface if we move along the orbits of and normal 
since we have gabW°‘^^ = 0 for any vector on the surface. 

If N corresponds to some (j) =constant surface, the corresponding normal one-form can be expressed 
as £q = Ada4>, where A is a scalar. Since we have assumed the four-dimensional metric gab to be 
non-degenerate, A 7 ^ 0. As proved in Appendix A.3.1, we will then have 

rVa4 := K = rdailnA) - y5,(r4) (13) 
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Note that while we assume = 0 on the null surface, it can in general be non-zero off the surface. So 
K has two contributions: one arises because A ^ 1 and the other because ^ 0 oS the surface. Very 
often, we will set A := 1 on the null surface making £a '■= a4> on it- The above result shows that k can 

be still non-zero on the null-surface due to the second term in its expression in Eq. (13). 

Since a£h '■= i^£b is a geodesic equation we see that the orbits of are geodesics. We may call k as 
the non-affinity coefficient because k = 0 if we choose an affine parametrization. Since lies on the null 
surface, the null surface can be thought of as being filled by a congruence of null geodesics. (If e.g., £“ is 
the Killing field normal to the horizon of a stationary black hole in an asymptotically flat spacetime, we 
can identify k with the surface gravity of the horizon [20]). Two specihc parametrizations of the metric 
on a general null surface are provided in Appendix C. 

3.2 Structure of the boundary term 

To study the boundary term on a null surface, it is more convenient to use the notion of a surface gradient 
rather than normal to the surface and rewrite the variation of the Einstein-Hilbert action in the following 
form: 

16^<5.4eh = - f d^x^G^^6gab+ f d^xv^^ - g^’^STl^)] , (14) 

Jv JdV 

We have again converted the total divergence VaW°‘ in Eq. (8) to a term on the boundary and defined Vc 
to be the surface gradient (not the unit normal) to the surface that bounds the volume. If the surface is 
given by the equation (j) = constant, then Eq. (14) is valid in any coordinate system where (j) is one of the 
coordinates with Vc = '£s^dc4>, the appropriate sign chosen as per the conventions of Gauss’ theorem (see 
Appendix A.l). When the volume integral is done as usual with the higher value of (j) at the upper limit 
of the integration, dc4> or —dccf is used depending on whether cj) or —cj) increases on going from inside 
the volume to outside through the surface. (Note that ^/—g appears in Eq. (14) since we have used the 
surface gradient instead of the unit normal. This is appropriate at this stage since V could be null in 
which case there is no natural notion of ” unit” normal ric or If choose to use the unit normal 

Ua, then na^/\^ will replace Vay/—g; see Appendix B.) 

We have earlier defined the useful function (5[Ac] in Eq. (9) for any one-form Ac as 

Q[Ac] = Acig^^^STl, - g^^SV:,) = V„((5A1) - (5(2V„A“) + (V^A,) Sg^^^ (15) 

where (5A“ = (5A“ -I- g°‘^SAi,. Adding a ^/—g and manipulating, we can obtain 

V^QlAc] = V^VclSAl] - 2(5(v^VaA“) + ^/^(VaAb - gcbVcA^)5g^^ , (16) 

When Ac = Vc = dc4> is the surface gradient, the integral of this expression gives the surface term of 
the Einstein-Hilbert action on a </> = constant surface, irrespective of whether the surface is timelike, 
spacelike or null. The surface term y/—gQ depends only on the surface gradient, da4>, at the 4> = 4>q 
surface. It does not depend on the normalization we may choose for Vc, i.e, it does not depend on B if 
we choose to write Vc = Bdc4>- H also does not depend on the behaviour of Vc away from the surface. 
But the decomposition in Eq. (16) needs to be rewritten if we choose to work with Vc = Bdccj). Also, the 
decomposition cares about the behaviour of Vc away from the surface due to the presence of derivatives 
of Vc- 

Let us now concentrate on the null-surface and denote the surface gradient on it by Vc = (-c = c4>- 

Following the route we took for non-null surfaces, we shall first convert the first term in Eq. (16) into 
a surface term on the 3—surface (f> = (f>o. For this, we require the notion of surface covariant derivative 
on the null surface, which, in turn, requires a projector on to the surface. Based on the analogy with 
non-null surfaces (where hab = gab — enaUb is the projector), if we try hab = Qab — it does not 

work because now hab£°‘ = £b — e(-°‘£a£b = 0, since £a£°‘ = 0 on the null surface. Thus, there is no 

straightforward extension of the projector in the non-null case to the null case, essentially because the 
metric is degenerate. 

The solution to this problem was suggested by Carter in [32]. Since we are not able to define a 
projector by using la alone, we choose an auxiliary vector such that iak‘^ = —1 at (j) = ‘Po- Then, as 
we prove in Appendix A.3.2, pab = gab + lakb + Ibka is the object on the null surface analogous to hab 
(see also [7]). Since is under our control, we shall assume that lak°‘ = —I and = 0 to be valid 
everywhere. Thus, we have 


lat := 0; = -1; = 0 . 


(17) 
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We shall enforce the first relation only on the boundary null surface. Further, we shall demand that 
these relations are respected by the variations. Since la = da4>i this demand on the first relation requires 
^g4><t> — This is a restriction on the variations of the metric and translates to the statement that we 
are only considering variations of the metric that keep the null surface null. We should ideally consider 
arbitrary variations of the metric, but the geometrical setup that we are using will lose its validity for 
the varied metric if the boundary surface ceases to be null. We hope to address this issue in a future 
publication. The second and third relations do not require any constraint on the metric variations. This 
is because these relations are the only restrictions on and we can choose an appropriate k'^ = k'^ + Sk°‘ 
in the varied metric. (Compare this with the case in usual double null formulations where /c“ is also 
chosen to be hypersurface-orthogonal [33].) Our analysis in Section 3.3 will be general while we shall 
choose a form for ka in Section 3.4. 

3.3 Boundary contribution in terms of 2-surface variables 

In order to find an expression for the boundary term on a null surface in terms of the 2-surface variables, 
similar to the expression for non-null surfaces in Eq. (12), we shall start with the expression for the 
surface term in Eq. (16), with Va = da4> being the surface gradient. We may take our null normal to be 
la = Ada4> = where A is some scalar. Then, recalling the definition of Q in Eq. (9), we can write 
the boundary term on the null surface in the following form: 

V^QlVa] = ^Q[la] = ^ {V^ValSll] - 26{^Wan + ah - gat^cn^g'^^} ■ (18) 

But since we do not seem to have any natural way of fixing the factor A (unlike the case in the non-null 
case where the condition riaU® = e was a natural choice) we shall make the choice A = 1. (The counter¬ 
term for a general A is derived in Appendix G.) This offers two advantages in our manipulations. First, 
it eliminates the 1/A factor in Eq. (18). Secondly, we will have Sla = daScj) = 0 as the scalar (f) labelling 
the null surface is not being varied. With the choice A = 1, the boundary term on the null surface 
becomes 

V^Q[h] = V^^c[Sil] - 2Si^Van + ^TgiVah - gabVct)5g’^^ . (19) 

Our first task will be to separate out a surface term from the first term in Eq. (19). We shall label this 
term as We have 

V^Ql = [Sil] = da [v^«l] (20) 

Let us consider this expression in the coordinates (</>, j/i, 2 / 2 , 2 / 3 ), where ( 2 /ij 2/2j 2 / 3 ) s-re some coordinates 
introduced on the =constant null surface under consideration. Then, the partial derivatives with 
respect to ( 2 / 1 , 2 / 2 , 2 / 3 ) s-re on the null surface and the (/-derivative is off the surface. Thus, we have 

V^Qi = d4^si{] + d^y/^sin ( 21 ) 

where a runs over ( 2 / 1 , 2 / 2 , 2 / 3 )- Now, S£^ = 6£°i_la = S {la£°')- We have assumed this to be zero on the 

null surface but not off it. Hence, the first term in Eq. (21) is not zero and ^/—gQl is not a total surface 

derivative on the null surface. 

In order to separate out a total surface derivative on the null surface from ^/—gQl^ there are two 
projectors, and q^, that we can use on (5£“ . They are defined by (see Appendix A.3.2) 

+ khb , ( 22 ) 

= 6^ + k‘^h + l^kb . (23) 

Now, n“j, projects orthogonal to h but projects orthogonal to both £a and ka- Since our aim is to 
separate out derivatives along the surface (i.e, orthogonal to h) even if they are not orthogonal to ka, 
we shall use the projector n“j,. Thus, 

V^Qi = ^a[^^^6£l] = da[V^U\5£i] - 

= da[V^Il\5£i] - y/^5 {lain [(5 ( 44 )] 

= da[V^U\S£’i] - ^k^db [5 ( 44 )], ( 24 ) 

where the last step was obtained by using our assumption S (lah) = 0 on the null surface. The first 
term in Eq. (24) is a surface derivative on the null surface as = 0. The second term in Eq. (24) is 



a bit of a bother. With (.a = dafj), S (£o^“) = Hence, this term has variations of the derivatives of 


the metric. We shall take out the 6 to obtain 

- k^db [<5 (4r)] = [k^db (4r)] + sk^db {lat) . (25) 

The vector 6k°' lies on the null surface as iaSk°‘ = 6 {iak°‘) = 0 as per our assumption 5ia = 0. Since 
= 0 on the null surface, the last term in Eq. (25) is zero and we are left with 

-k^db[5{ian] = -5[k'‘db{tan] . (26) 

Substituting in Eq. (24), we obtain 

V^gQi = da[^gn\5t^] - [k^db (4r)] 

= da[V^ii\se\] - s [V^k^db (£ar)] - ^ [k^db (4r)] g,, 6 g^^ (27) 

Thus, we have written yJ—gQi in a form where all the variations of the derivatives of the metric have 
been smuggled into surface derivatives and total variations. 

Substituting Eq. (27) back in Eq. (19), we find the following expression occurring in two places: 

Vat + ^da + kHbVd^ = ■ (28) 

Thus, the boundary term on the null surface reduces to 

V^Q[ia] = da [V^U%6ii] - 2S(^U%Vat) + ah " gab^^ct)5g^^ . (29) 


We have succeeded in separating out a total derivative on the surface and a total variation to remove all 
derivatives of the metric. In order to understand what degrees of freedom need to be fixed on the null 
surface, we shall now decompose into < 7 “^, h and 

First, note that the relations q°^^£a = 0 and q‘^^ka = 0 are respected by the variations. This is because 
we have assumed Eq. (17) to be valid even on variation. Thus, terms of the form £ahSq°'^, £akbSq°'^ and 
kakbSq°'^ would reduce to zero. For example, 

£akbSq'^^ = S (g“''4) kb - q‘^'^kb8£a = 0 . (30) 

Using this result, we can simplify gab8g‘^^ as follows: 

9ab6g^^ = gab [Sq^^ - S (tk^) - 6 (tk^)] = qab6q^^ + 2 {£akb + £bka) 5 (tk^) 

= qabSq‘^^ + 2£bka6 (tk^) = qabSq‘^^ - 2kaSt, (31) 

where we have also used £a5k°‘ = 0 , arising from our assumptions of 5£a = 0 and 5{£ak°‘) = 0 . 

Next, we shall simplify {Va£b) Sg‘^^. This is a symmetric object since £b = Vb4>- We have 

{Va£b) 5g^'^ = (ya£b) 5q^^ - 28 (rfc") 

= (Va4) - 28£‘^k^Va£b - 28k^£‘^Va£b 

= (Va£b) 8q^^ - 28£^kWa£b - 2n8k^£b 

= {Va£b)8q‘^^ -28£‘^k^Va£b, (32) 

where we have used £°‘Va£b = K£b (see Eq. (13)) and £a8k°‘ = —k°‘8£a = 0. The expression {Va£b) 8q^^ 
can be simplified as follows: 

(Va4) < 59 “' = {Vm£u) 8q^^ = {qT - £^ka - k^£a) {q^ - £^kb - k^£b) {Vm£n) 8q‘^^ 

= {qT - £^ka) (g," - Tfcb) {VM ^9“' 

= {q^q^Vm£n - qatkbVm£n " qb£^kaVm£n) 8q°-^ 

= {Qab - 2Kq^kb£m) 8q^^ 

= Qab8q^'^ (33) 
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where we have used ia5q°'^ = <5 [q°‘^ia) = 0 to get to the second line, kak^Sq^-^ = 0 (see Eq. (30))to get 
to the third line, the definitions of 0a& and k (see Appendix A.3.5 and Eq. (13) in this paper) and the 
symmetry of Va(-b to get to the fourth line and £ 3 ( 7 “^ = 0 to get the final result. 

Substituting back in Eq. (32), we get 

(Va4) < 55 “^ = 0a6<5<7“'’ - 2 {k^yjb) 5t . (34) 

We shall also write n“^Va£^ in the following form: 

= {qt - rfcfe) Vat = 0 + K, (35) 

making use of the definitions of 11“;, and k and introducing the expansion scalar 0 (see Appendix A.3.5). 

We have thus derived all the results that we need. Substituting Eq. (31), Eq. (34) and Eq. (35) in 
Eq. (29), we get 

V^Q[L] =da - 26 [v^(0 + k)]+^ [0a& - (0 + l^)qab] 6q‘^^ 

+ 2^^ [(0 + K)kc + i‘‘Vckb] dt (36) 

Thus, the surface term in the action for a null surface takes the following form: 


SAnull = / d^X y/^Q[la] 

Jav 

= [ d^x {\/^Va (n“t,(5£')_) - 26 [v^(0 + k)] + [©ah - (0 + n)qab] 6q' 
Jav 

+2^[(0 + Ac)fc, + £^Vc/cb]<5r} . 


(37) 


In the first line, we have a structure like we have in the non-null case in Eq. (164) (in Appendix B), 
reproduced below for convenience: 




non—null 



^/\h\Da{6nl) + 6{2./\h\K) - VW\{Kab - KKb) 6h^^ 


(38) 


The first term in Eq. (37) is the total 3-derivative term (since 11'^^ = n“j,£a = 0), the second term can 
be cancelled by a counter-term and the third term is the term that may be killed by fixing q°‘^ on the 
surface. We have a crucial extra term with 61^^ which we shall deal with when we discuss degrees of 
freedom. 

The result in Eq. (37) is invariant under coordinate transformations on the three surface. That 
is, it is invariant under coordinate transformation from (^,to ((/),x'^,x'^,x'^). The integra¬ 
tion volume i/—gd^x is invariant once (j) is fixed since yf—gd^xdcj) is invariant under general coordinate 
transformations, while the integrand is a scalar once £a = Va<(> and a particular is chosen. 

The result in Eq. (37), with its ^/—g and d^x, may seem unusual for readers with some familiarity to 
the literature on null surfaces. To remedy this, we shall write this result in a specific coordinate system 
((/), A, adapted to the null surface. This will bring the integration element to the form given, 

e.g., in [7] (also see Appendix A.3.3). In the process, we shall also choose a form for fc^. The double 
null coordinate system commonly found in the literature (see, e.g., [33]) will be a special case. This 
coordinate system, as well as the meaning of our constraints on variations in this coordinate system, is 
described in the next section. 


3.4 Writing the Results in A Special Coordinate System 

In this section, we shall detail the construction of a coordinate system adapted to a null surface and also 
choose a form for ka- This coordinate system has the double null coordinate system [33] as a special case. 
In fact, the form of the metric on the null surface will be identical to that in double null coordinates. 
Please note that we require the form of the metric only on the fiducial null surface and that is what we 
shall write down. The expression for the metric in this section is valid only at 0 = and should not be 
taken to represent the form of the metric in a region around the null surface. 

The object qab has been constructed to represent the induced metric on a 2-dimensional space orthog¬ 
onal to £a and ka- Let us choose /c“ such that ka := VaA = 9aA for some A at the boundary null surface 
(p = (pQ. This is always possible. (The GNC coordinates, whose construction is given in Appendix C.I, 
is an example. In these coordinates, £a = daX and ka = —daU.) Like (p, A is not varied when we vary 
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the metric and hence 5ka = da{5\) = 0. This is an additional restriction which was not present in the 
derivation leading up to Eq. (37). Note that the relation ka = 9qA is assumed to hold only at ^ = (/iq. 
If it is assumed to hold everywhere, then we are specializing to a double null foliation. But that is 
not required for our purpose. (Again, the GNC metric is a good example, ka = {2ra, —l,rj3^) in the 
coordinate system (u, r, x^) and is equal to —daU on the null surface r = 0.) In this set-up, qat represents 
the induced metric on the 2-surfaces at the intersection of our fiducial (j) = (j)Q null surface and the many 
A =constant surfaces (whose normal coincides with ka at <f> = (jjo). If we choose coordinates (^ 1 , 2 : 2 ) on 
a particular 2-surface and carry these coordinates along integral curves of £“ on the null surface, then 
the inverse metric on the null surface in the coordinate system ((( 1 , A, 2 ^, 2 ^), using = 0, kak°‘ = 0, 
= —1, valid on the null surface, and the fact that both and are orthogonal to the basis vectors 
on the (j) = (fo, A =constant 2-surfaces (which means that the metric is block-diagonal), can be written 


(39) 


down as 


( 

-1 

0 

0 



^ab 

-1 

0 

0 

0 



9 = 

0 

0 

gii 





1 0 

0 

gi2 

9^^ 

/ 

giving the metric. 

on the null surface (f = 

(fo, in 

the form 






-1 

0 

0 

\ 



-1 

0 

0 

0 



gab = 

0 

0 

911 

9l2 




^ 0 

0 

gi2 

922 

/ 


(40) 


where qn, < 712 , etc. have the appropriate expressions in terms of etc. For this metric, the 

determinant satisfies —g = g, where g is the determinant of the 2-metric qab- The relation between 
the determinant of a 2-metric to the determinant of the 4-metric in a general case has been derived in 
Eq. (65) in Appendix A.2 as 

9 = g^^gW _ (gA0)2 ’ ('^1) 

where we have inserted A instead of the 1 that stood for the coordinate in the appendix. Then, since 
we are keeping g°^^£a£b = g^^ = 0, g°‘^kakb = g^^ = 0 and g°^^iakb = g'^^ = —1 fixed even under variation, 
Eq. (41) gives 

\/^ =^/q, <5 ^ (v^) • (42) 

Note that our constraints fix g^^ and g^^ while the other components are free to vary. In particular, 
these variations may disturb the relations £°‘ = —dx°‘/dX and k°‘ = —dx°'/d(f) valid for the on-shell metric. 
In these coordinates, we can write the boundary term in Eq. (37) as 

5AnuU = / dX(fz± y/q Q[ea] 

Jav 

= [ dXd^Z^ {^Va - 26 [^{Q + k)] +y^[Qab " (O + K)qab] 6q^^ 

Jav 

+2Vg [(0 + k ) fcc + £^^ckb] 5t] . (43) 

Here, d'^z±_ = dz^dz^. 


3.5 Analyzing The Results 

The final decomposition of the boundary term in the variation of the Einstein-Hilbert action is given in 
Eq. (37) and it is written in specific coordinates and a particular choice of ka in Eq. (43). We see that 
the action principle has again suggested what can be done! On a null surface, we can add the integral 
of 2 y^^(0 -I- k) {2ytq[Q + k] with the coordinates and choice of ka of Section 3.4) as the counter-term. 
Further, we need to set 5q°‘^ = 0 and = 0 on the boundary, which requires fixing the set (g“^,^“) (six 
independent components, since = 0) on the boundary. This is analogous to fixing (six independent 
components) on the non-null surface. To know for certain whether this furnishes a well-posed variational 
problem, we need to analyse the structure of the equations of motion on the null surface. As we have 
mentioned in the introduction, there is no standard way of writing the Einstein equations on null surfaces. 
Hence, we postpone this analysis for a future work. We shall briefly discuss about the degrees of freedom 
in Section 4. 
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The expression Eq. (43) also tells us that the momentum conjugate to g®*” is ^/q [©at — (0 + K)gab] 
and the momentum conjugate to is 2y/q [(0 + k) fee + . We can write the momentum conjugate 

to in another form as 


= -51{Q + k)\ . (44) 

In the cases where k = 0, as would happen if = 0 everywhere, the counter-term would be the integral 
of just 2y/qQ. In this case, we can make use of 0 = {1/y/q){d^Jq/dX) (see Eq. (123) in Appendix A.3.5) 
to write the counter-term as 


2 f dXd^z± y/qQ = 2 / dXd^z±_ ^ d^z±y/q 

JdV JdV JddV Ai 


2[S{X2)-S{X,)] , 


(45) 


where Ai and A 2 are the limits of integration of A. Thus, the counter-term in this case, apart from an 
overall factor, is the difference in the area of the 2-surface orthogonal to £ and k on the null surface, 
denoted by 5'(A), at the limits of A integration. A similar interpretation can be given to the counter-term 
for the non-null case, the integral of —2^J\h\K over the boundary surface, if we assume that the Lie 
bracket of the normal vector n“ with any of the coordinate basis vectors on the surface is zero (see 
Appendix E). This, in turn, requires to be a tangent vector to a geodesic congruence, which may not 
be true in general. With this assumption, —y^\h\K is the rate of change along n“ of the volume of a 
3-surface element orthogonal to n“. 

We have provided a decomposition of the boundary term on the null surface working purely on the null 
surface. This decomposition can also be obtained by treating the null surface as a limit of a sequence of 
non-null surfaces. In the next section, we shall demonstrate that a particular method of taking the null 
limit of Eq. (46) gives back the decomposition in Eq. (37). 


3.6 Null surface as a limit of sequence of non-null surfaces 

The previous sections provided a decomposition of the boundary term on the null surface working purely 
on the null surface. In this section, we shall show that this decomposition can also be obtained by 
treating the null surface as a limit of a sequence of non-null surfaces. 

Let the relevant null surface be part of a family of surfaces, not all of which are necessarily null, 
characterized by the constant values of a scalar function The null surface is specified by (j^ix) = tpo- 

(That is, the surface (f) = is null but the other (p =constant surfaces may be spacelike, timelike or 
even null.) We shall label the null surface as M and also introduce the scalar function r = p — pQ. Thus, 
we can represent the null surface by the symbol Af or by one of the two equations, (() = ^0 or r = 0, as 
per convenience. The normal to the surfaces has the form £a = AVap, for some scalar function A{x). 
If the surfaces were spacelike or timelike, we could have fixed A{x) by demanding = ±1. But now, 
we have with = 0 at p = pQ. Thus, imposing this constraint would make A diverge 

at p = pQ. Hence, we cannot impose any such normalization condition without introducing infinities. 
The expression for the normalized normal for a non-null surface is given in Eq. (141) in Appendix B: 
Ub = N'S/bp- On the other hand, we have used the expression la = dap in obtaining Eq. (37). So, 
we have the identification Ua = Nla- Now, N is defined in terms of in Eq. (142) in Appendix B: 
g’l’l = e/N"^ . So, the null limit can be imposed by demanding that N —>■ 00 . Near the null surface, we 
shall take g^^ to have the behaviour = Hr" -|- 0(r"+^) for a finite non-zero B and some positive 
integer n. 

We shall now examine the decomposition of the surface term of the Einstein-Hilbert action for a 
non-null surface, provided in Eq. (164) in Appendix B, and obtain the null surface limit of each of the 
terms. For convenience, we reproduce Eq. (164) below: 




dV inon—null 



^\Dai6nl) + S{2^\K) - ,/\h\iKab - Khab) 


(46) 


In obtaining the limit, we shall assume that the metric components are finite and that the metric 
determinant is finite (and non-zero) when we take the null limit. (Later on, we will demonstrate these 
conditions and the final result in specific coordinate systems.) 

Let us start with the surface boundary term in Eq. (46). This term has been analyzed in Appendix F.l 
to show that it can be written in the form: 


^\Da {5nl) = da, [A/=^<5r] + 2da, [^l^^SiaN] . 


(47) 
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The first term in Eq. (49) is finite in the N ^ oo limit. In the second term in Eq. (49), we use the 
expansion for near the null surface and see that the r factor cancels (as we have taken the variation 
to not affect the coordinates) and the term is finite in the null limit of iV —?► oo or r —>■ 0. Thus, we have 
decomposed \/\h\Da (^n“ ) in such a way that each term in the decomposition is clearly finite in the null 
limit. 

On the other hand, the surface boundary term for a null surface in Eq. (37) has been shown in 
Appendix F.l to reduce to 

da = do. [V^sr] (48) 

Thus, from Eq. (47) and Eq. (48), we obtain the result that the null surface limit of the surface boundary 
term in Eq. (46) is 


{5nl) da [y/^iv\5£\] + 2do [^r<51niV] . (49) 

Next, let us look at the counter-term in Eq. (46). From Appendix F.2, a similar analysis shows that 
the null limit of the variation of the counter-term turns out to be 

5{2^J\h\KY=°-25[^{Q +K)\-25[^f^ i°‘do\riN] . (50) 

Consider now the last set of terms in Eq. (46). Among them, the hijSK^^ term is manipulated in 
Appendix F.3 leading to 

Vlft) ’'=° - (0 -b k) QabSq'^^ + 2yf^ (0 -b At) kbSt - 2yf^ (0 -b At) 5 In TV 

- In N) g.jSg'-^ - 2^g (r In A^) 5 In A^ (51) 


Again, one can see that each term is finite in the null limit under our assumptions. 

The last term to be considered in Eq. (46) is the term. The pertinent expression is derived 

in Appendix F.4 as 


- ^/\h\K,j5h^^ = + 2Sefy,kj + 2 {d. In N) St + 2e {d. In N) S In N 


(52) 


Adding Eq. (49), Eq. (50), Eq. (51) and Eq. (52), we obtain the following result when we take the null 
limit of the integrand on the right hand side of Eq. (46): 


r—>-0 


y/\h\Da{Snl) + S{2y/\h\K) - VWliKab - Khab) Sh^’^ 
do - 25 [v^ (0 + At)] + [0aA, - (0 + At) gab] 

-b 2y/—g5t [(0 -b At) /ci -b ^^Vikj\ 

-b 2do In A^] - 25 [yf^ rdo In A^] - 2y/^ (0 -b At) 5 In A^ 


- In N) g,jdg^^ - 2^ {rdo In A^) 5 In A^ 


yf^ 2 {di In N) St + 2t {di In A^) 5 In N 


(53) 


We find that the first two lines reproduce the result in Eq. (37). Thus, we should be able to show that 
the rest of the terms cancel out. This is proved in Appendix F.5 and we obtain the following expression 
for the surface term: 


[5A, 


1 —^■O 


dV inon—null 


= [5A 


av\null 


f d^X {^Va [n“,5^i] - 25 [yTg (0 + «)] + ^ [0ah " (0 + «) qab] Sq^’’ 

Jav 

+2y/^St [(0 -b At) A:, + tV^kj] } . (54) 


where we have replaced the surface index a in the first term in Eq. (53) with the four-dimensional index 
a as = 0. Thus, we have re-derived the result in Eq. (37), with proper canonical momenta conjugate 
to reaffirming our faith in the correctness of this result. 

In Appendix C, we have illustrated the result in Eq. (37) or Eq. (54) for two specific parametrizations 
of a general null surface. 
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4 Discussion of the Number of Degrees of Freedom 

We shall now try to enumerate and identify the metric degrees of freedom that we have to fix on 
the null boundary in our approach. The number of degrees of freedom that we have to fix on the 
boundary modulo gauge invariance (i.e, after making use of gauge invariance to eliminate certain degrees 
of freedom) corresponds to the number of physical degrees of freedom in the theory for a well-posed action 
(see Section 2 and [2]). The addition of a boundary term to the action can change the number of degrees 
of freedom to be fixed on the surface. But, for a well-posed action, the number of degrees of freedom to be 
fixed on the boundary should match with the number of degrees of freedom obtained from the analysis of 
the initial value problem of the equations of motion. For an initial value problem, the number of degrees 
of freedom is half of the number of initial data to be specified. For a variational problem, the number of 
degrees of freedom is the number of boundary data to be fixed at the boundary. It is well-known from 
the analysis of the initial value problem of general relativity, with an initial spacelike surface, that the 
theory has two physical degrees of freedom per spacetime point (see, e.g., [19, 20]). The same conclusion 
can be reached by analysing the variational problem for a non-null boundary. In the next paragraph, we 
shall demonstrate this for a spacelike boundary. The analysis for a timelike boundary is similar. 

Consider the case when the boundary is a spacelike surface in a 3 -I- 1 framework. Initially, the metric 
has ten degrees of freedom but we have the freedom of making four coordinate choices. We shall effect 
a 3 -I- 1 split by slicing the spacetime into t = constant surfaces with the fiducial boundary being a 
t = constant surface. The t-coordinate of each spacetime point is fixed. Next, we introduce the normal 
to the t-constant surfaces. There are four components to the normal Ua- But since it is the normal to 
t = constant surfaces, Ua = —Ndat and there is only one functional degree of freedom in Uq in the choice 
of the function N. Imposing the normalization condition na'n°‘ = — I fixes N in terms of the metric 
by the relation = — 1/7V^. This relation is preserved on variation since we take the varied normal, 
Ua + 6na, also to be normalized. The variation of the normal is Sua = —{SN)dat, since coordinates are 
kept constant during the variation, and essentially contains the variation 6g**. On the other hand, the 
upper components are n“ = —Ng‘^’’dat = —Ng°‘^ = {\/N,—Ng°^*), where we have substituted for g** in 
terms of N. From this relation, we observe that the time evolution vector = (1,0, 0,0) can be written 
as -I- A^“, for the purely spatial shift vector with components N^g‘^* (for more details, 

see Chapter 12 in [I]). Thus, n“ = (1/iV, —iV“/iV) and i5n“ basically contains Sg** and the three Sg*°‘. 
Introducing the notation for the spatial components of the metric gap, the degrees of freedom in 
the 3-I-I decomposition can then be said to be in (TV, N°‘,hap) or equivalently in (n“, hap). That is, the 
ten degrees of freedom in the metric have been inherited by the normal (four degrees of freedom) and 
the induced metric hap (six degrees of freedom). Out of these ten, variations of only six, the degrees of 
freedom in hap = gap, appear in the boundary term in the variation of the Einstein-Hilbert action, after 
we remove the counter-term part and the surface term on the 3—surface, and hence only these need to 
be fixed on the t = constant boundary surface. But we have three more coordinate choices to make. 
The choice of three coordinates on the initial t = constant surface can be used to remove three of the 
six degrees of freedom in hap on the surface. This leaves 3 degrees of freedom. It is generally argued 
that one of these three degrees of freedom just gives the location of the 3-surface in the 4-dimensional 
spacetime, as determined by the slicing of the spacetime into t =constant surfaces, and can be tuned by 
appropriate choice of the t coordinate [19, 20]. Thus, only two degrees out of the ten degrees of freedom 
in the metric are real physical degrees of freedom and only these need to be fixed on the boundary if we 
use our freedom of choice of coordinates appropriately. 

Coming to the case of null surfaces, the initial value problem on a null surface has been analysed by 
several authors [22, 23, 24, 25, 26, 27, 28, 29, 30, 31]. But unlike the case for the initial value problem 
on a spacelike surface, there appears to be no one standard formulation for null surfaces. The constraint 
structure and the initial data to be specified varies with the formalism. Thus, the analysis of the question 
of whether the action is well-posed with the counter-term that we have proposed is not as straightforward 
as in the case of non-null surfaces. We intend to pursue this question in a future publication. 

For the time being, we can ask which degrees of freedom in the metric have to be fixed to set our 
null boundary term in Eq. (37) to zero. But even this question cannot be satisfactorily answered in our 
framework. The degrees of freedom to be fixed on the boundary appear to be in the three components 
of (coefficient of has only three components since it is symmetric and orthogonal to and 
fc“) and the three components of = iaSi"" = 6 = 0 as per our assumptions). This gives 

six components and one might think that four components might be eliminated by diffeomorphisms to 
match with the two degrees of freedom that one obtains in a non-null case. But notice that we have 
imposed the constraints £°'£a = 0, = 0 and £ak°' = —1 throughout, even under variations. The 
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effect of the last two constraints on the variations of the metric can be specified only once is specified. 
But, with ^a = dafj), S = 0 imposes the constraint Sg'£'£ = 0. In Section 3.4, we take ka = daX on 

the null boundary to obtain the other constraints as = 0 and Sg^^ = 0. Hence, in our framework, 
variations of certain degrees of freedom of the metric are constrained. To do a rigorous analysis, either 
the framework has to be extended to eliminate these constraints or these constraints have to be explicitly 
taken into account (by adding Lagrange multipliers, for example). 


5 Conclusions 

Our aim in this work was to find out what boundary term can be added to the Einstein-Hilbert action 
for a null boundary to make the action well-posed. We wanted to do this from first principles, adhering 
to the philosophy: The action can itself suggest the counter-term to be added to make it well-defined 
and also what boundary conditions go along with this counter-term. In the case of the Einstein-Hilbert 
action with a null boundary, we expect the action itself to suggest what counter-term needs to be added 
to make the action well-defined. 

We first undertook a general analysis of the surface term for a null surface (p = (po, or r = (p — (po = 0, 
with la = dap and being the auxiliary null vector satisfying k^'ia = — 1. We took p as one of the 
coordinates and considered variations that vary the metric but do not affect the coordinates. (This 
assumption also means 5£a = 0.) The variations were taken to respect the following constraints: i) t.t := 
0, ii) fe.fe = 0 and hi) (..k = — 1 , with the first constraint being imposed only on the null surface. In 
this case, we found that the surface term on the null surface can be expressed in the following form (see 
Eq. (37)): 

SAnuU = I d^X {- 2S + k)] + [0a6 - (0 + K)qat] Sq‘^^ + PcS£^} . 

JdV 

(55) 

where, taking the normal £a = dap for the null boundary being & p = constant surface and to be the 
auxiliary null vector, we have = S^-\-k°-£i,, $£']_ = 5£°‘ g°‘^5lb and qab = gab + £akb + ka£b, the induced 

metric on the null surface. Qab = qlffqb'^m£n is the second fundamental form, 0 = 0 “ is the expansion 
scalar, k is the non-affinity coefficient on the null surface and Pc = 2%J—g [(0 -\- k) kc-\- £^'X/ckb] ■ (The 
definitions of Qab, 0 and k can be found in Appendix A.3.5 and Appendix A.3.1.) In a coordinate system 
adapted to the null surface and with a specific choice of ka, we can write the above decomposition as 
(see Eq. (43)) 

6Anuu = I dXd^z^ {^Va {Tl\5£i) - 26 MQ + k)] + ^ [0^6 - (0 + K)qab] 6q^'^ + Pc6£^) , 

JdV 

(56) 

The integration is over a parameter A, varying along the null geodesics such that ka = VqA locally, 
and two coordinates (z^,z^) that are constant along the null geodesics, q is the determinant of the 2 - 
metric qAB- The conjugate momentum to turns out to be ^yq [0q;, — (0 -I- K)qab] while the conjugate 
momentum to £“ is Pc = 2ytq [(0 q- k) + £^'^ckb] ■ In Appendix G, we have repeated the whole analysis 
for a null normal of the general form £a = Adap to obtain a decomposition with the same structure as 
in Eq. (55) or Eq. (56). We have also verified the result in Eq. (55) using two parametrizations of the 
metric near an arbitrary null surface, the GNC and NSE metrics, in Appendix G. 

We note that the structure is very similar to the structure in the non-null boundary (see Eq. (46)). 
The first term in Eq. (55) or Eq. (56) is the total 3-derivative term (due to the properties of H^^), the 
second term can be cancelled by a counter-term, 2 ^/^{Q n) or 2 y/q{Q n) respectively, added to the 
action and the remaining terms can be killed by fixing q°‘^ and £°' on the surface. 

The question of whether this prescription leads to an action that is well-posed requires answering 
the question whether the suggested boundary conditions and the equations of motion are consistent and 
whether there is a unique solution for any choice of the boundary conditions. This task is complicated 
by the plethora of frameworks available for dealing with Einstein equations on a null surface. Even if 
the action turns out to be well-posed, the analysis of degrees of freedom of the theory on the null surface 
requires careful analysis of the constraints that we have imposed. There are six degrees of freedom in 
( 7 “^ and which matches the number of degrees of freedom found in hab for a non-null boundary. But, 
unlike the standard analysis for non-null boundary, our framework constrains certain variations of the 
metric. For example, the constraint that the null surface remains a null surface means Sg^^ = 0. There 


15 


is also the matter of whether we can find a better method to normalize the normal on the null surface 
than the A = 1 fiat that we imposed in this paper. Detailed investigations in these directions have been 
kept for a future work. 

We summarise in Tab. 1 the similarities and differences between non-null surfaces and null surfaces 
with respect to the intrinsic geometry and treatment of the boundary term. 


Table 1: Comparison Bet-ween Non-null and Null Surfaces. 


Properties 

Non-null 

Null 

Normal 

Ua, nan°- = e = ±1 

o 

II 

“Dimension” 

3 

2 

Induced Metric 

13 

qAB = gabe^eB 

Auxiliary Vector 

-en“ 

= -1 

Integration Measure 

d^Xy/\h\ 

d^xdXy/q 

Second Fundamental Form 

Kab = -h'^Vmnb 

Qab = qTQb'^mrin 

Counter-term 

-/d^xvW 2K 

f dXd^ZxyAj 2 (0 -I- k) 

Boundary Term 

^^\Da5n°^ 

yhjVa (n%wi) 

Degrees of Freedom 


qAB, 

Conjugate Momentum to Induced Metric 

-^/\h\ 

[0ab _ qab (0 

Conjugate Momentum to Normal 

0 

-v^ kb - 26^, (Va^)] 


The defining characteristic that distinguished a non-null surface and a null surface is the surface 
gradient, whose norm is non-zero for a non-null surface and zero for a null surface. The norm of the 
surface gradient for a non-null surface can be normalized to ±1 by multiplying with a suitable scalar. 
The induced metric on a non-null surface is a 3-metric ha/s while the induced metric on a null surface is 
a 2-metric qab- Hence, the “dimension” of the null surface has been listed as 2. For a null surface, we 
define an auxiliary vector to be such that = — 1 and k°'ka = 0. For a non-null surface, the term 
auxiliary vector is not generally used but we have put —en“ in the respective place in the table as this 
vector satisfies (—en“) Ua = —1. The rest of the entries in the table are quantities which we have already 
discussed in the analysis of our results on the decomposition of the boundary term on a null surface. 
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Appendices 

A Some Requisite Pedagogical and Background Material 

A.l Gauss’ Theorem 

In order to explain the conventions regarding the use of the Gauss’ theorem, we refer to the proof of the 
theorem given in Chapter 3 of [7]. The Gauss’ theorem is stated in the following form: 

[ d^X^VaA^^ = [ dJ:aA\ (57) 

Jv JdV 

with dTja being the directed surface element on the boundary of the integration volume. Note that [7] 
uses Greek indices to run over all spacetime indices and Latin indices to run over indices on a surface, 
opposite of the convention we have adopted. Then, the proof is illustrated using a set of =constant 
surfaces. The following expression is obtained on integrating over a volume between = 0 and = I 
surfaces: 


d'^Xy/^VaA' 


d^xV^A° 


(58) 










This may be rewritten as 




d^3 




i:° = l 




V^gA^ 


r:°=0 


= ® d^Xy/^A°-(.a 


(59) 


with la = at a:° = 1 and ia = —daX^ at = 0. This depends on the fact that the integration 
is carried out from a;® = 0 to a;® = 1. If the integration was done the other way, the signs would have 
flipped. 

So we may formulate the following rule. Let us assume that x^ is going to be integrated from lower 
value to higher value. In that case, dax'^ or —daX^ is to be used at an =constant surface according to 
whether x^ or —x^ is increasing as we move from inside the integration volume to outside through the 
surface. 


A.2 Decomposition of y/—g in Terms of the Determinant of Metric on a 
2-surface 


One relevant question is whether there is a decomposition of \/—g in terms of y^, q being the determinant 
of the 2 -metric qab on the null surface, akin to the decomposition \/—g = N^/\h\ in the timelike 
and spacelike case. We shall prove in this appendix a preliminary result on the decomposition of the 
determinant of a 4 x 4 metric in terms of the determinant of a 2 x 2 submatrix, which will be later applied 
to a null surface. 

We start by writing down a general result relating the determinant of a 2 x 2 submatrix to the 
determinant of the whole matrix. We shall prove the result working with the metric written in the 
coordinates x^, x^, x^) with the components on the (p =constant surface being denoted by hap (refer 
Appendix B): 


/ g4>4> g^^ g<i>3 \ 

501 ^11 hi2 hi3 

502 ^21 ^22 ^23 

\ 503 ^31 ^32 ^33 / 


(60) 


In this case, we can use the definition of an inverse matrix element applied to to write 



(61) 


where h is the determinant of hap, the 3 x 3-matrix obtained by deleting the ())-column and 0-row from 
gab- Now, we can play the same game again with hap. The determinant of the 2x2 matrix qAB, 
A,B = 2, 3, defined by qAB = hAB satisfies an analogue of Eq. (61): 


h = 


g 

/ill ’ 


(62) 


where is the 11-th component of the matrix /i“^, the inverse of the matrix hap. Substituting for h 
in Eq. (61), we obtain 


^ g 
9 g'P'Ph'^i 


(63) 


Now, the denominator above can be expanded as follows: 




(64) 


which is easiest to obtain by using the formula /i“^ = 5 “^ — en°‘nP. Thus, we obtain a relation relating 
the determinant of a 2 x 2 submatrix with the determinant of the full 4x4 matrix: 


^_g_ 

3 g 00 gll _ (gl 0)2 • 


(65) 


A.3 Null Surfaces 

In this section, we shall discuss various pedagogical and otherwise useful material on null surfaces that 
we will be needing in the main text. 
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A.3.1 The Non-affinity Coefficient k 


In this appendix, we shall prove that i°‘Vah oc £b, for £a being the null normal, £a = Adacj), to a null 
surface (j) = (j)o. We have 


rVa4 = £‘^Va{Adb^) = t^-^daA + rAVaVf,^ 

= r + £<^AVbVa4> = t\daA + t A\l bi!^^) 

= t^-^daA + rVf,4 = [£‘^da{lnA)]£b + ^dbi£‘^£a) (66) 

Consider some coordinate system with tp as one of the coordinates, say (</>, cci, a; 2 , 2 : 3 ). In such a coordinate 
system, 

£a=Adap=iA,0,0,0) . (67) 

Now, db{£°‘£a) will only have the <))—component at the null surface. This is because £“£q = 0 all along 
the null surface and hence only d^{£°'£a) 0 at the null surface. Thus, we obtain db{£°‘£a) oc £b and 

rVa4 = ( 68 ) 

where k is a scalar. It may be termed the non-affinity coefficient as it will be zero for an affine parametriza- 
tion of the null geodesics [34]. An explicit expression for k is derived in Appendix A.3. 6 . 

A.3.2 Induced Metric on a Null Surface 

In this appendix, we shall discuss how to find the induced metric on a null surface. If we choose an 
auxiliary vector such that k°'£a = — 1 on the null surface (f> = po (note that we have not specified k°‘ka 
yet), then we have the following results on the null surface. First, consider the object 

n% = 51 + k<^£b . (69) 

It is easy to verify that = 0, so any vector can be acted on by 11^^ to give = II^^jL^ such 

that AI°'£a = 0. To verify that an operator P is a projector, we need to verify = P, which is satisfied 
in this case as = II'^^j. 

While may be good as a projector, hab had the additional status of being the induced metric 
on the surface. Now, flat, does have the property that if we look at the components on the surface 
(represented by Greek letters, meant to run over in a coordinate system {cj),x^,x‘^,x^)), they 

do satisfy Hab = g^/s as £a = 0. But if we want a symmetric object, we may turn to 

Qab — 9 ab £akb “b k(i£b — Il^^b “t“ £akb • (^ 6 ) 

This also has the property q^£a = 0 and, in addition, q^£^ = 0 at ^ (while £l‘l£^ = £“). To see if 
this is a projector, consider 

9b 9c = 9c + ^“fc&9c • (71) 

Thus, in order for q^ to be a projector, we need the additional condition £°‘kbq’p. = 0, or kbq^ = 0. This 
means 

fcb9c = kc- kc + kbk’’£c = 0 , (72) 

which requires k°'ka = 0 , i.e, we need to choose to be a null vector. 

Thus, for a null surface with normal £a, we choose an auxiliary vector fc“ such that (i) £ak°‘ = —1 
and (ii) = 0 . Then, the symmetric object qab is such that q^ acts as a projector to the space 
orthogonal to In fact, it also acts a projector to the space orthogonal to fc“. Thus, it projects to a 
2-dimensional subspace of the tangent space on the 3-surface. We may say that it is a projector to the 
2-surface orthogonal to both and fc“. The 3-dimensional space orthogonal to £°‘ a,i (p = po corresponds 
to the tangent space of the null surface under consideration while a further projection orthogonal to k°‘ 
takes us down to a 2-surface on the null surface. To formalize this statement, let us define a set of 4 basis 
vectors (£“, e‘\), with A = 1 , 2 , at each point on the null surface. e‘\ are a pair of linearly independent 

spacelike vectors which satisfy £ae‘X = 0 and kae‘\ = 0. We shall call such a basis as a canonical null 
basis, and introduce the notation 

V(a) = (Afe,e^), (73) 
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with, e.g, V( 3 ) = ei. Now, we have 


= 0; = 0; q^e\ = e‘X ■ (74) 

Thus, q^ projects to the subspace spanned by e'\. The metric induced on the surface formed by the e'\ 
is 

^akb ^bka) — Qah^A^B — QAB • (75) 

Now, qAB contains the whole information about qab since ga6^“ and qabk^ are already constrained to be 
zero. Thus, qAB, and hence qab, represents the induced metric on the two-dimensional surface spanned 
by ca and orthogonal to £ and k. 

To get the induced metric on the 3-surface, note that the vector space on the 3-surface is spanned by 
(£,6^) as £“£q = 0, e'\£a = 0 but /c“£q = —1. Thus, the induced metric on the 3-surface consists of the 
components = qAB, gab£‘^£^ = 0 and gab£°‘e% = 0. Hence, the 3-metric is also effectively qAB- 

Let us now introduce the dual basis [20] to the canonical null basis. We need a set of four linearly 
independent one-forms, vl“\ such that Denoting the inverse of the 2-metric qAB by 

we introduce a set of two one-forms, e^, such that 

ef = g^^gabCs = q^^{eB)a ■ (76) 

Multiplying both sides by qAC and we get the inverse relation 

. (77) 

Then, it can be easily checked that 

vl“) = (-k,-^, e^), (78) 

with k representing the one-form with components ka etc., provides the required dual basis. In particular, 
we have 


e^ef = e\q^^{ec)a = q^'^gabec^A = g^^QCA = ^a ■ (79) 

The canonical null basis and the dual basis allows us to write down the following decomposition of the 
Kronecker tensor: 

51 = -tkb - fc“4 + e\et , (80) 

which can be easily checked by contracting the lower index with each member of the canonical null basis 
and the upper index with each member of the dual one-form basis. Raising the lower index, we obtain 

= -tk^ - k^-t + 

= -tk^ -k^£^ + q^^e\e%, (81) 


which implies 

Contracting both sides with and e^, we can invert the above relation to obtain 


ab _ „AB a b 

g — g ^A^B 


r^CD 


ab„C„D 


g = 9 e, 


(82) 

(83) 


If we lower one index in Eq. (80), we can use Eq. (77) to obtain 

gab = gABeaef , (84) 

the inverse of which, qAB = qab^A^’s’ gives us back our original definition of qAB hr Eq. (75). 
So, finally, let us list out all the relations between the 4-dimensional tensors, qab and 5 “**, and the 
2 -dimensional tensors, qAB and q^^: 


qAB = qabSA^B , 
A B 

qab = qABSa Sf, , 


,AB 


g 


— g ^a^b , 
AB a b 

=g e^eB ■ 


(85) 
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A.3.3 Erecting a Coordinate System on the Null Surface 

Suppose we have a set of coordinates x°' = charting the four-dimensional spacetime with 

the null surface under consideration. Any set of three continuous, infinitely differentiable functions, say 
?/^), of the spacetime coordinates x“ constitutes a system of coordinates on the null surface 
provided the set of values of these functions at each and every point on the null surface is unique. Then, 
the coordinate basis is the set of three vectors 


a dx^ 

“ ■ 


( 86 ) 


If gab be the components of the metric of the ambient spacetime in the coordinates x“, the induced 
metric on the null surface is given by [7] 

hap = gabe^ep ■ (87) 


This is a 3-metric and the determinant of this 3-metric, h, will be zero as the surface is a null surface. 
The easiest way to see this and to work on the null surface is to erect a coordinate system naturally 
suited to the null nature of the surface [7]. The null surface is filled by a congruence of null geodesics, the 
integral curves of the normal vector We choose a parameter A varying smoothly on the null generators 
such that the displacements along the generators are of the form = £“dA. (In this paper, the null 
surface is taken as a (j) = constant surface and the normal is taken as ia = Ada4> for some A. £“ is then 
fixed and A has to be chosen appropriately. If our aim is just to erect a coordinate system on a given 
null surface, we may first choose A to be some parameter which varies along the null geodesics and then 
choose £“ = dx“/dA.) If we further ensure that A varies smoothly for displacements across geodesics, we 
may choose it as one of the coordinates on the null surface. The other two coordinates are to be chosen 
as two smooth functions 2 ^ = ( 2 :^, 2 ;^) that are constant on each null geodesic. They act as a unique 
label for each null geodesic. In this coordinate system, varying A would correspond to a displacement 
along a particular null geodesic while varying the set would correspond to displacements across the 
generators along points of equal A. The basis vectors in the coordinate system (A, z^) are 


e 


a 

\ 


dx°‘ 


= 


i9x“ 


A = 1,2 . 


( 88 ) 


Note that the identification = dx°'/dX is possible because the coordinates have been chosen to be 
constant along the null geodesics. Previously, we had chosen an auxiliary vector and then demanded 
that e\ satisfy e‘\ka = 0. But while we are working purely on the null surface, there is no notion of an 
auxiliary null vector, as no vector on the null surface will satisfy k°'ia = —I- In fact, having defined the 
coordinate basis vectors on the null surface, if we now we want to move into the ambient 4-diinensional 
spacetime and define qab, etc., we can specify uniquely by the four conditions k°‘ka = 0, k'^ia = ~1 
and gabk°^e^A — *-*• (This is how we defined /c“ in the case of GNC metric in Appendix C.l.) 

The components of the induced metric are h\\ = gab^’^i’^ = 0 (since £“ is null), h\A = gabf^°‘e\ = 0 
(since e\ lies on the surface and is the normal) and qab = gabeAS^B- coordinate order (A, z^, z^), 

the matrix form is 

/ 0 0 0 \ 

hab ~ j 1 ■ (89) 

\ 0 qi2 q22 ) 

The determinant of hab, h, is obviously zero and will remain zero even if we transform to a different 
coordinate system on the null surface as the determinant changes only by a Jacobian factor. The line 
element is 2-dimensional: 

= qABdz^dz^ . (90) 

Keeping this coordinate system as a reference, we can explore other coordinate systems on the null 
surface. The metric in any coordinate system = {y^, , y^) on the surface is given by 


,, dx'^ dx^ 

- qAB gyp 


(91) 


In general, none of the components need to be zero, although the determinant will vanish. But consider 
the special case (?/^(A, z^,z^), 2 /^(z^,z^),?/^(z^,z^)). This coordinate system has coordinates y^ and y^ 
constant on the null geodesics, y^ may now be considered as the parameter varying along the null 
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geodesics and we will have = M{y^,y^,y^)dx°‘/dy^ (with M = dy^/dX). In this case, the metric will 
again take the form 

/ 0 0 0 \ 

Kb=[ 0 gii q[2 ■ (92) 

\ 0 q'i 2 922 / 


A.3.4 Directed Surface Element for the Null Surface 

For any set of coordinates y°‘ = {y^,y‘^,y^) on the null surface with associated basis vectors e“ = 
dx°'/dy°‘, the invariant directed surface element for the null surface is given by [7] 

dEa = eabcde\ele'^dy^dy‘^dy^ . (93) 

Now, £ 060^616263 is proportional to the normal ia (for a null surface, there is no unique normalized 
normal, but our treatment here will work for any choice of the normal) as its contraction with any vector 
on the null surface, expressible as a linear combination of e“, is zero. Thus, we may write 

£ 060^616263 = /4 (94) 

for some scalar function /. Contracting with we obtain 

/ = -eabcdk°‘e\ele’^ . (95) 

Since / is a scalar, it can be evaluated in any coordinate system. So we choose the 4-dimensional 
coordinates to be (<)>, y^,y‘^, ?/^) with (j) constant on the null surface. The form of the normal will then 
be (.a = Ada4> for some scalar function A. In this coordinate system, borrowing the notation = from [7] 
to indicate equalities valid in a specified coordinate system, 

/ = -£0123fc‘^ = ^ , (96) 


where we have used the condition k°'ia = k^iff, = — 1 to find k^. If we choose the eabcd tensor so that 
£0123 is positive, then we shall have 


f ^ \/^ 

J ^ , 

with the choice ia = da4> used in our treatment of the boundary term leading to 

f = V^9- 


(97) 


(98) 


Thus, for the null surface represented as cj) = constant and the normal being specified as ia = Ada4>, th® 
directed null surface element, in a coordinate system where (j) is one of the coordinates, becomes 

dEa = ^^^^iady^dy'^dy^ = y/^da(l)dy^dy'^dy^ . (99) 

Note that, while this result is not generally covariant, it is valid in any coordinate system with (j) as one 
of the coordinates. In the case of ia = da4>, it reduces to simply 

dEa = V^iady^dy'^dy^ . ( 100 ) 

In fact, nowhere in the derivation of Eq. (99) and Eq. (100) did we have to make use of the assumption 
of the surface being null. In the case of non-null surfaces, we have the extra luxury of specifying the 
normalization of ia uniquely to give Ua = ia = Ndacj) (see Eq. (141)). The auxiliary vector should 
then be taken as k°' = —erA as the only condition we demand of in the above derivation is k°'ia = — 1 .. 
Then, A = N and ^—gjA = and we arrive back at the familiar result for surface element for a 

non-null case. 

We have mentioned that Eq. (99) and Eq. (100) are valid only when ip is taken as a coordinate system. 
One way to formulate a fully covariant expression for the surface element is as follows. Note that all 
the treatment till Eq. (95) is fully covariant. The non-covariance sneaked in at the evaluation of /. 
Now, given ia = da(p and the triad of vectors e“, we can formulate a tetrad of vectors (—fc“,e“), with 
= (—1,0,0, 0) in (^, 2 /^,?/^,?/^) coordinate system. Thus, satisfies = —1 but is not necessarily 
null. Using this tetrad and the metric, we can construct ten independent scalars = gab{—k°‘){—k^): 
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g<pi = 9 ab{—k°')e\, etc. In the coordinate system (0, , these scalars are just the components of 

gab- But now imagine taking the matrix corresponding to the metric gab and replacing each component 
with the corresponding scalar to get a matrix of scalars. The determinant of this matrix is also a scalar. 
Let us denote it by g. Hence, from Eq. (98), we have 

f = ^ f = ^ f = (101) 

which is an equality between scalars and hence valid in all coordinate systems if valid in one. Note that 
this also means 

y/^ = -eabcdk^'elelei , ( 102 ) 

from Eq. (95). 

In the special coordinate system (A, z^) introduced in Appendix A.3.3 with ia = da4>, the directed 

surface element reduces to [7] 


dT,a = ia\/qdXdz^dz^ , (103) 

where q is the determinant of the 2-metric qab- Note the difference in minus sign from the analogous 
expression in [7]. This minus sign difference arose because [7] defines the normal as £a = —da4> so that 
1 °" is future-directed. 

A.3.5 Second Fundamental Form Qab and Expansion Scalar 0 

In this section, we shall introduce the second fundamental form Qab and the expansion scalar 0 for a 
null surface. The terminology follows [34]. But [34] works with a foliation of null surfaces while we have 
only specihed the (j) = (f>o surface to be a null surface among the (j) = constant surfaces. Hence, we cannot 
blindly carry forward the results in [34]. But we shall see that the results that we require are in fact 
unaltered. 

Following [34], we use the projector in Eq. (69) to introduce the extension of the second fundamental 
form for the null surface at any point, with respect to the chosen normal £ (since, unlike the non-null 
case, the normalization of the normal is not fixed and hence the choice of £ is not unique), to the tangent 
space of the four-dimensional spacetime manifold at that point: 

Qab = . (104) 

Now, consider the following object on the null surface: 

dadt^cid = (ffa + ^"ka) {U\ + £'^h) Vdd 

= Qab + + edaU^VJd + tkai^Vcld 

= &ab + n\dhVdd + fikaU\£d + Kka£‘^kb£d 

= Oab + lli\kbVc {£^£d) , (105) 

where we have made use of the results £d°‘b£a = 0 and £“£q = 0 on the null surface. The second term 
in Eq. (105) contains an expression of the form Vb{£°'£a)- Let us simplify this expression to a form that 
will be more useful to us. We write 

Vb(r4) = gbcVdtla) 

= idbc - ibkc - kb£c)Vd£da), (106) 

where we have used Eq. (70) in the second line. In Eq. (106), the third term has the combination 

rVe(r4) = £^ddt£a) = 0, (107) 

since ddc is a derivative along the null surface and the value of dla is zero throughout the null surface. 
Next, consider the first term in Eq. (106), which has 

gj,cV=(r4) =qldd£da) 

=qtd^{tta), (108) 
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in the coordinate system in which (j) is one of the coordinates, since the only direction along which £°'ia 

varies on the null surface is off the surface i.e the direction along which (j) varies. 

Now, using Eq. (67), 

qt cx q^lb = 0 . (109) 

Thus, qf is identically zero on the null surface and the first term in Eq. (106) also vanishes. (To prove 
that a tensor quantity is identically zero, it is sufficient to prove that it is zero in one coordinate system.) 
So, we arrive at the result that 

V6(r4) = -rve(r4)4 . (110) 

Using this result in Eq. (105), we obtain 

0a6 = q^qi^cid ■ (111) 

Instead of labouring forth with ” the four-dimensional extension of the second fundamental form”, we 
shall henceforth take the liberty of referring to Qab as just the second fundamental form of the null 
surface. To be precise, we should also add that it is the second fundamental form with respect to the 
normal ia under consideration but we shall consider this to be understood henceforth. 

The Qab is a symmetric object. To see this, we use Eq. (67) and decompose Vah as 

Vaib = AVaVb(l) + ebda\nA . (112) 

The first term in Eq. (112) is symmetric while the second term does not contribute to Qab- Hence, 
proved. 

We shall next prove a relation between Qab and the Lie derivative of qab along £. The Lie derivative 
formula is 

£iqab = dqab + qad^b^'^ + +qdb^a£‘^ ■ ( 113 ) 

Substituting for qab from Eq. (70), we have 

£iqab = {£a^ dkb + kb^ d£ a + dka £ ka^ d^b) b^a + {£akd£ ka£d)^ b£^ a£b£ {£bkd£ k^id)^ a£^ ■ 

(114) 

Contracting with q'^q^, and using q^ia = 0 and q^ka = 0, we obtain 

q^ql£eqab = q^Qn aib + Vh^a) = Qmn + Qnm = 20mn (115) 

which implies 

Qmn = ^qm<l^£t<lab ■ (116) 

Finally, let us look at the trace of Qab- 

0 = 5“'0ab = g“'’0ab . (117) 

We shall refer to 0 as the expansion scalar on the null surface. More specifically, it is the expansion 
along of the 2-surface on the null surface orthogonal to The reason for this terminology becomes 
clear if we take the trace of Eq. (116). We obtain 

0 = ^q‘^^£iqab - (118) 

This equation may be further manipulated to a form easier to interpret. We can make use of the canonical 
null basis introduced in Appendix A.3.2. We have 


qab = QABeaSb , 

ab _ „AB a b 

q — q ^A^B - 


where summation over A, H = 1, 2 is implied. Then, 

0 = ^q^^e\e%£i {qcneaeb) = ^q^^ £eqAB + q^^ qBceA£iea 
= ^q^^£eqAB + eA£iea = ]^q^^£iqAB - ea£ee\ - 


(119) 


( 120 ) 
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where we have used e^ef = 5^, q^^qsc = <^c i^A^a) = "^^2 = 0 along the way. Now, 

let be the coordinate vectors e“^ and 6“^ in the special coordinate system (A, introduced in 

Appendix A.3.3. Since and are then members of a coordinate basis, their Lie bracket should be 
zero, i.e, 

[£,eA]=0^il^e^ = 0. (121) 

Enforcing this condition, we obtain 

e = ^ . (122) 

2 V? 


where q is the determinant of the 2-metric qAB- Thus, we see that 0 represents the fractional change in an 
area element of the A = constant 2-surface on the null surface as it is Lie dragged along We can replace 
the Lie derivative in the above expression with an ordinary derivative since qAB = qab^A^^B ^ scalar 
under 4-dimensional coordinate transformations with the kept as the same physical vectors. Thus, 
£iqAB = (■°'daqAB which becomes £iqAB = dqAB/dX for the parameter A introduced in Appendix A.3.3. 
Hence, 


1 dy^ 

y/q dX 


(123) 


A.3.6 K and 0 in terms of (£, k) 

In this appendix, we shall derive expressions for the non-affinity coefficient n and 0 in terms of £a and 
auxiliary vector ka, for ia being the null normal, ia = Ada(f>, to a surface cj) = 4>o- 

Let us first define the auxiliary vector (need not be null for manipulations in this section) such 
that k°‘ia = — 1. Then, 


= rVaiAdb^) = tjdaA + tAVa^bCi^ 

= t^-jdaA + tAVbya<P = t^-jdaA + 1AV/-A) 

= t^-^daA + tVda = [tda{lnA)\£b + \db{tia) (124) 

Substituting in Eq. (124), we obtain 

tVJb = [rda(ZnA)]4 - ^d,{£Ha)£b 

= [£<^da{lnA) - ^d,{l^la)]£b (125) 

Hence, we see that a£b = i^£b with the non-affinity coefficient k given by the formula 

K = r5a(/nA) - ^d,{t£a) (126) 

Using this expression, we can also evaluate 

- rfc''Va4 - tk'^yatb (127) 

= 0 + K - . (128) 

Let us define 

R = -^da{£Hb) = -k'^tyjb , (129) 

a measure of how much iai^" varies as we move away from the null surface. If Vaff; is symmetric, like 
in the case when we have it = Vf,^, we shall get k = k. If Vo^6 is antisymmetric, i.e if ib is a Killing 
vector, then we get k = —k. Eq. (128) now takes the form 

Var = 0-bK-tk. (130) 


24 




A.3.7 Decomposition of y/—g in Terms of ^/q for a Null Surface 


We shall use the result Eq. (65) which is applicable even when we take the null limit. In the limit the 
(j) = constant surface under consideration is null, g'^'^ = 0 and h = 0. Taking the limit g’^’^ —>■ 0 on 
Eq. (65), we get 


a = 


-q 

(gl0)2 ■ 


(131) 


For the null surface, we need q to be the determinant of the 2-metric qAB = gabe^e’^g (see Eq. (75)). 
To apply result Eq. (131), we shall specialize to a coordinate system such that e\ are coordinate basis 
vectors and la = da4>- On the null surface, we shall introduce the two coordinates , constant 

on the null geodesics and a third coordinate g, which is a parameter varying along the null geodesics 
such that the vectors e'\ lie on a g = constant surface and 6(4 = dx°‘fdz^. We shall also have = 
(1/M) dx°‘/dg for some scalar M. In the case where we choose the special coordinate system (A, 
introduced in Appendix A.3.3, we will have M = 1 and = dx°‘/dX. In the coordinates ((/, /r, , 

gAB = qAB, g^A = gab (r/M) = 0 and g^f, oc gahl’^l^ = 0. Also, consider gabl°'k^ = iaxb/M) Now, 

gabl°'kk’ = iak°' = k’^. Hence, only contributes to gabl°'k^- Since gabl°'kk’ = —1, we have k'^ = —1 and 
we obtain g\b = M. Thus, in the coordinates [(j), g, z^, z"^), the metric takes the form 


^ 500 

M 

501 

502 

\ 

M 

0 

0 

0 


501 

0 

511 

5l2 


V 502 

0 

5l2 

522 

/ 
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Using Eq. (131) (or by direct calculation of the determinant from the above matrix), we can write down 


a=-M\ , 


and 

\M\y/q . 

If we specialize to /i = A, we would have g = q and 

\/^= y/q ■ 


(133) 

(134) 

(135) 


B Boundary Conditions for Spacelike and Timelike Surfaces: 
An Alternate Approach 

In this appendix, we shall detail an alternate approach to arrive at the standard prescription for hxing 
the boundary conditions on timelike and spacelike surfaces [18]. As we have seen in Eq. (8), the surface 
term of the Einstein-Hilbert action is the integral of the quantity 

yTgQivc] = V^Vc(g‘^’’Sr:, - g^>^5T:,) (136) 

over the boundary of the spacetime region under consideration. Here, Vc is the surface gradient. To 
be concrete, let us take <f) = (po to be the surface, with pQ a constant. Then Vc = idc4>, with the sign 
decided according to the conventions of the Gauss’ theorem (see Appendix A.l). Let us assume that 
our surface is such that p increases on going from inside the integration volume to outside through the 
surface. Hence, we shall use Vc = dcp from now on. It is clear that Eq. (136) contains the variations of 
the metric and its derivatives. We can separate out the terms with the variations of the metric and its 
derivatives in Eq. (136) and write 

V^Q[Va] = Va - ^V^P^‘^‘^^5{ddgab)] , (137) 

where = 1/2 {g^'^Sa — g^^^Sa). Since Vc = dcp, the coefficient of 5{ddgab) is —We can 
easily check that ^ 0 in general. Hence, Eq. (136) contains variations of the normal derivatives of 

the metric. 

Now, let us consider the case of timelike and spacelike surfaces as the boundary and see how we can 
fix boundary conditions without having to fix both the metric and its normal derivatives at the boundary. 
Generalizing Eq. (136), we shall define Q[Ac] for any vector Ac to be 

Q[Ac] = Acig^^^STl, - g^^^SV:,) . (138) 
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For timelike or spacelike boundaries, we can normalize the surface gradient to obtain the unit normal Ua 
such that 

; non“ = e, (139) 

with e = 1 for a timelike surface and e = — 1 for a spacelike surface. We shall also demand that the 
normalization be preserved under the variation of the metric, i.e, 

5{nan°‘) = na5n°' + 5nan°‘ = 0 . (140) 

We choose a coordinate system with (j) taken as one of the coordinates. Let the other coordinates, the 
coordinates charting the (j) = constant surface, be labelled , x^). As mentioned in the introduction, 

the Greek indices (other than (p) run over {x^,x‘^,x^) while the Latin indices run over (cj), x^,x'^,x^). Let 
us use N for the normalization factor in Ua, i.e 

ria = Nda4> = Nva , (141) 

where N is to be assumed positive so that Ua and da4> point in the same direction. The normalization 
nan-“ = e then relates N to the metric as 




g<k<k = g/7v2 ^ (^42) 

Another convention that is used often is to demand that p increases in the direction of the normal vector 
n“, i.e., n°^dap > 0. In this case, we would write 


fia = eNva, (143) 

with N positive, where we have put the tilde just to distinguish the normal in this convention from the 
Ua in Eq. (141). Now, from the standard procedure of calculating the inverse of a metric we know that 

^00 ^ Cofactor [g 00 ] ^ Pet [fee./;] ^ h 

where the 3x3 matrix corresponding to the metric on the p = po surface, obtained by deleting the 
p-TOw and (((-column from the matrix representing the metric, was denoted by ha/s and the symbol h was 
defined to be the determinant of this matrix. We can rewrite Eq. (144), using Eq. (142), as 

g = eN'^h = -N'^\h\ . (145) 

Thus, we arrive at the following decomposition for y/—g'. 

= N./\h\ . (146) 

Using this expression and Eq. (141), we can rewrite the expression for the surface term in Eq. (136), in 
the case of timelike or spacelike surfaces, as 

V^Qivc] = VW\nc{g‘^'’Sri, - 

= ^A^Q[nc] ■ (147) 

As we can see from Eq. (137), this expression contains the variations of the metric as well as its normal 
derivatives. Our aim is to discover a counter-term such that, when added to the Einstein-Hilbert action, 
the surface term obtained in the variation will contain only the variations of the metric. That is, we 
would like to express the surface term in the form 

= 5[A] + YabSg'^^ . (148) 

Then, it is —X that we have to add to the action as the counter-term. 

Proceeding towards this goal, let us first manipulate (3[Ac] for an arbitrary vector A‘^. We have the 
following relations: 


S{XaAb) = XaSAb - Sri,A, ; (5(VaA“) = XJA^ + . (149) 
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Using these relations, we can rewrite the expression for Q[Ac] in Eq. (138) as follows: 

= Va{SAl) - (5(2VaA“) + WaAb 8g-\ (150) 

where 

5A1 = . (151) 

In the hrst line of Eq. (150), we have used the notation [A, B] for the commutator, AB — BA, and {A, B} 
for the anticommutator, AB + BA. 

Let us now specialize to the case of Ua- Then, we have 

Q[nc] = Va(^nl) - 5(2Van^) + V5g^\ (152) 

where (5n“ = SrA + g°‘^6nb. Using Eq. (140), we can see that 

Snlua = 0, (153) 

which means that the vector 5n“ lies on the surface (p = 4 >q. This property can be used to decompose 
the first term of Eq. (152) in the following manner: 

Va{5nl) = 5tVa{5n\) 

= {hi + en°-nb)S/ a{5n\) 

= hlVa{8n\) + eif-UbV a{8n\) 

= Da{5nl_) - e(n“Vanb)5n5_ 

= Da{5n°[_) — eab8n\ . (154) 

On the way to obtaining the above expression, we have introduced the induced metric 

hab = gab - eUaTlb , (155) 

and a derivative operator Da such that DaV’’ = For vectors U“ on the surface (i.e satisfying 

V’^Ua = 0), Da is the natural covariant derivative on the 3—surface (see [7] and Chapter 12 in [1]) 
compatible with the induced surface metric. 

To proceed further, let us look at the nature of the variations we are considering. The 8s are field 
variations on the metric and leave the coordinates unchanged, which is why they could be taken inside 
integrals over spacetime and derivatives with respect to the coordinates. The scalar (/) is taken to be 
fixed during the variation, which means that the foliating surfaces (j) = constant are kept fixed during 
the variation. In our case, this is not an extra assumption as we have already taken (j) to be one of our 
coordinates. Now, from Eq. (I4I), 

Ua = Nda4> Sria = 8Nda<i) = 8{\a.N)na . (156) 

This shows that 8na is in the direction of Uq. We can write this in terms of variations of the metric using 
the constraint rAua = e, equivalent to Eq. (142). We see that 

,51niV = ^8\nN^ = -ijlng'^'^ = -^8g‘^‘^ = -^8g‘^^nanb . (157) 

Thus, we obtain the following explicit expression for the variation of na in terms of the metric: 

8na = -^8g'‘^ninjna ■ (158) 

Now, the second term in Eq. (154) contains 

ab8nl_ = ab8n^ + a^8nb = abna8g°'^ + ab8nag°'^ = abna8g°'^, (159) 

where we have first used the definition of 8n\ from Eq. (151), and then applied the result a!’8nb = 0, 
which we know from Eq. (156), twice. Thus, Eq. (154) can be rewritten as 

Va((5n“ ) = Da{8n\) - eabna8g°'^ . (160) 
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So, Eq. (152) takes on the form 


QH = DaiSnl) - (5(2V„ri“) + (VaUb - eUaUb) Sg‘^^ (161) 

Our manipulations have naturally produced the quantity VaTib — snaUb- This quantity is nothing but 
the negative of the extrinsic curvature [1] defined by 

Kah — c^b — ^a^b “t“ ^’bldCLb i 

where we have used the definition of the induced metric given in Eq. (155). Kab has the following 
properties: 

K,, = Kj,- =0; K = K,, = -V^n^ (163) 

The second equation tells us that KijSg^^ = KijSK^^. Using this result and the results in the last three 
equations, we can rewrite the surface term of the Einstein-Hilbert action, in the following useful form: 


[5Agv] 


non—null 



\/\h\Q[nc] 



^/\h\Da{Snl) + S{2./\h\K) - 2KS{^/\h\) - \/\h\Kab 5h‘^^ 

^/\h\Da{Snl) + S{2^\K) - ^\{Kab - Khab) ■ 


(164) 


On integrating over the 3—surface, the first term in Eq. (164) is a 3—divergence which can be converted 
to a boundary term. If the surface bounding the integration volume is a closed surface, then this term 
goes to zero. The second term in Eq. (164) can be killed off by adding the counter-term —2^J\h\K to 
the Einstein-Hilbert Lagrangian. The last term in Eq. (164) can be put to zero by fixing the induced 
metric on the boundary. 

If we choose to work in the convention where the normal is chosen to be the one in Eq. (143), then 
we would have the surface term in the form 


i^A^vLon-null = f f^VW\Q[nc] 

JdV 


dV 

= f dPx 
JdV 


iy/\h\Da{5h\) + 5{2ey/\h\K) - ey/\h\{Kab - Khab) Sh' 


ah 


(165) 


C Boundary Term on The Null Surface for Specific Parametriza- 
tions 

In this appendix, we shall verify the validity of the result in Eq. (37) by working out two specific cases. 
We shall introduce two parametrizations for a general null surface and calculate the boundary term of 
the Einstein-Hilbert action on the null surface. In both cases, we find results in agreement with Eq. (37). 


C.l Gaussian Null Coordinates (GNC) 

The first parametrization that we shall discuss is what is commonly known as Gaussian null coordinates 
(GNC), in analogy with the Gaussian normal coordinates [19]. Gaussian normal coordinates are con¬ 
structed by extending the coordinates on a non-null hypersurface to a spacetime neigbourhood using 
geodesics normal to the surface. This prescription does not work for a null surface as the “normal” 
geodesics lie on the null surface itself. Hence, the construction of Gaussian null coordinates is carried out 
by making use of certain, uniquely defined, auxiliary null geodesics. As far as we know, this coordinate 
chart in the neighbourhood of a null surface was introduced by Moncrief and Isenberg in [35] and hence 
may also be termed as Moncreif-Isenberg coordinates. The construction of GNG is also described in [36], 
[37] and [38]. 

C.1.1 Construction of the Coordinates 

We shall now describe the construction of GNC. Our description will be less technical and more intuitive. 
We shall mostly follow the notation of [38]. Note that [38] has the GNC metric in Eq. 3.2.7 but then 
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claims that it can be further constrained and writes it as in Eq. 3.2.8, i.e makes the replacement a —>■ ra. 
We do not subscribe to this claim and will use Eq. 3.2.7. 

Consider a smooth null surface in a spacetime manifold M. Let gab represent the metric on M. 
We shall take the spacetime to be four-dimensional but the following construction can be easily extended 
to an n-dimensional Lorentzian spacetime as long as the extra dimensions are spacelike. Take a spacelike 
2-surface, C, on N and introduce coordinates on that surface. We have shown that the null 

surface N is generated by null geodesics. The null geodesics cannot be along the spacelike surface at 
any point, since any vector that lies on the spacelike surface has to be spacelike. Thus, one necessarily 
moves away from the spacelike surface by travelling along any of these null geodesic. Therefore, the null 
geodesics can be used to define coordinates on the null surface as follows. Let u be a parameter, not 
necessarily affine, along the null geodesics with u = 0 on ^ and increasing towards future. Let each null 
geodesic be labelled by the coordinates {x^,x^) of the point at which it intersects C- Then, any point 
on in a neighbourhood of sufficiently small so that the geodesics do not cross, can be assigned the 
coordinates (u, x^), where (x^, x^) corresponds to the label given to the null geodesic passing through 

that point and u corresponds to the value of our chosen parameter at that point. Henceforth, when we 
talk about, say, coordinates on the null surface, it is to be understood that we are talking about such a 
neighbourhood. We shall call the future-directed vector field tangent to the null geodesics, d/du, as {. 
and the basis vectors corresponding to the coordinates (x^, x^) as = d/dx^ with H = 1, 2. So far, we 
have not introduced the fourth coordinate, and hence the above definitions, with the partial derivatives, 
will continue to hold only if the fourth coordinate is chosen to be constant on the null surfaces. But such 
a choice will indeed be made, as we see next. 

Having thus constructed a coordinate chart on the null surface, it is time for us to move out into the 
surrounding spacetime. We shall do this with the help of a new set of null geodesics. At each point on 
the null surface, there is a unique vector (note the interchange of notation of £ and k compared to 
[38]) satisfying the following four constraints: (i) kak°‘ = 0, i.e it is a null vector; (ii) £ak°‘ = — 1 and 
(iii) Xfjia = 0 for A = 1, 2. Erom £ak°‘ = —1, it is clear that this vector does not lie on the surface at any 
point but “sticks out”. Thus, we can make use of this vector to “go off” the null surface. From each point 
on the null surface, send out null geodesics in the direction of labelled by the coordinates, (u,x^,x^), 
of that point. The fourth coordinate may be chosen as an affine parameter along the null geodesic. Let r 
be that affine parameter such that r = 0 represents the null surface and k = —djdr. (The sole purpose 
of the minus sign is to reproduce the exact form of the metric in [38] by off-setting the fact that we 
demand iak°‘ to be —1 while [38] equates the same to -1-1.) These constraints uniquely determine the 
affine parameter. Any r' = A(u,x^,x^)r + B(u,x^,x^) is also an affine parameter, but the condition 
r' = 0 on the null surface will lead to H = 0 and the condition = — 1 for fc = —djdr' fixes the value 
of A to be unity. Having chosen this affine parameter, any point in a sufficiently small neighbourhood 
of the null surface (so that the null geodesics do not cross) can be assigned the coordinates (M,r, x^,x^), 
where the set (u,x^,x^) corresponds to the label of the null geodesic that passes through that point 
while r represents the affine parameter value at that point. 

Let us now turn our attention to the form of the metric in the neighbourhood of the null surface. We 
shall start by writing down the form of the metric on the null surface. Since i is null everywhere on the 
null surface, we can write down our first metric component: g^u = gab£‘^£^ = 0 , everywhere on the null 
surface. Next, since the null surface is represented by r = constant, = d/dx^ at the null surface has 
to lie on the null surface. This means that X\£a = 0, which implies gabX\£^ = guA = 0 everywhere on 
the null surface. Finally, denoting the rest of the components relevant for vectors on the null surface, 
gAB = gabX^Xg, by gAB, we have the following list of components of the metric on the null surface: 

guu = 0 ; guA = 0 ; gAB = Qab', (166) 


Let us next determine the metric components in the neighbourhood of the null surface that we are 
considering. First, since k = —djdr is null everywhere, we have gabk°‘kf' = grr = 0 everywhere. To find 
the rest of the components, let us write down the geodesic equation governing the null geodesics that 
are integral curves of fc. Since r is an affine parameter, we have 

fc“Vafc'' = 0 ^ = 0 ^ r,,, = 0, (i67) 


where we contracted with gab in the last step. Putting a = r tells us that grr remains zero along the 
geodesics. But if we put a = g, where x^ is any one of the three surface coordinates (u, x^, x^), we have 


L^rr — 2 ( A ‘^drgrfi) — drgrfi — 0 ■ 


(168) 
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Thus, grfi is also constant along the geodesics. (This is a common feature for affinely parametrized 
geodesics. For null geodesics, if d/d(j) is the tangent vector with (j) being the affine parameter, then we 
shall have drf,grf,a = 0. For non-null geodesics, we can get the same constraint provided (f) is taken as the 
proper time or proper length or a constant multiple thereof, so that we can put d^g^^ = 0 as in the last 
step above.) To fix the values of note that we have the constraints k°^^a = — 1 and kaX‘^ = 0 at 
r = 0, which translate into gm = 1 and grA = 0. From Eq. (168), we conclude that g^ = 1 and grA = 0 
everywhere in the region under consideration. Finally, defining qab = dabX'^X^, with the definition 
Xa = d/dx^ extended everywhere, and respecting Eq. (166), we have all the components of the metric 
as follows: 

grr — grA — 0; — 1; guu — 2rQ;, guA — ‘^I^Ai gAB — QAB : (169) 

where a and /3^ are finite at r = 0 so that g^u and gm vanish at r = 0. The 2x2 matrix qab is positive 
definite (since Xy 4 are spacelike vectors). The choice of the sign for g^A is not of importance as the sign 
can be flipped by changing the coordinates on the spatial surface, x'^ -A —x^. On the other hand, the 
sign of guu does have a physical significance. With this particular choice, the vector becomes timelike 
for r > 0 and spacelike for r < 0 when a > 0. Also, 

g’’’’ = 2ra + 

where fP = qabP^, which is a positive quantity since qab is positive definite. Thus, the normal to 
r-constant surfaces, daX = satisfies the condition 

r°‘ra = 2ra + r^/3^. 


The second term is strictly positive. But the first term dominates as r ^ 0. With a > 0, r-constant 
surfaces near r = 0 are time-like in the r > 0 region and space-like in the r < 0 region. In fact, all 
r-constant surfaces in the r > 0 region become time-like by this choice. 

So, finally, here is the line element in GNC: 

ds^ = —2radu^ + 2drdu — 2rPAdudx^ + qABdx^dx^ (II’O) 


Note that there are six independent functions (a. Pa, Qab) in this metric. We have used up the freedom 
of choosing the four coordinates to eliminate four out of the ten degrees of freedom in the metric. [The 
inverse metric and the Christoffel symbols corresponding to this metric are given in Appendix D.l. The 
Ricci tensor components are provided in [35].] Introducing a time coordinate by t = u -f r in place of u, 
we can rewrite this metric in the standard ADM form [21] as. 


ds^ = —N^dt^ — 2 (1 -I- ra) {dr — N^dt) {dr — N'^dt) 


2rPA {dr — N''dt) {dx^ — N^dt) + qab {dx^ — N^dt) {dx^ 
1 -I- 2ra + r'^P^ 

—2 — 2 ra + 




-rP^ 


3 -I- 4ra 


—2 — 2 ra + r 


2 02 


2 ra + 2 rPAN^N^ - 2(1 -f ra) (IV’')^ -f qABN^N^ 


N^dt) 


(171) 

(172) 

(173) 

(174) 


This expression shows that the relationship between the GNC metric components and the standard ADM 
variables {N,Na, ha/ 3 ) is not simple. In particular the degrees of freedom in /ia /3 comes in entangled in 
terms of the other degrees of freedom in GNC variables. 


C.1.2 Surface Term on a Null Surface as a limit in GNC 

Gaussian null coordinates (GNC) provide a situation where r = constant surfaces are time-like for r > 0 
and null for r = 0. In order to show the validity of our previous derived result for null surfaces, Eq. (37), 
we will obtain the boundary term on the null surface as the limit of the boundary term on a time-like 
surface, providing a specific example of the abstract manipulations carried out in Section 3.6. Consider 
the following expression for the term to be integrated on the boundary of a time-like or space-like surface 
as given in Eq. (46): 


= ^/\h\Da{Snl) + 6{2K^/\h\) - - Kh,,)5h^P (175) 
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where Q as given in Eq. (138) in Appendix B has to be evaluated for the unit normal. We shall use 
the object which, apart from a constant factor, represents the canonical 

momentum conjugate to in ADM formalism [1]. Then, the above equation will take the form 


Q^\ = VW\Da{5nl) + 5{2K^\) + (176) 

Further, we have the following relations 

P'^^hab = 2^\K- Pab5h<^^ = -P<^Hhab (177) 

Thus, we can also write the boundary term Q\/\h\ using the above expressions for P°‘^hab and Pab^h°‘^ 
leading to 

Q^/\h\ = y/\h\Da{5nl) + (178) 

The idea is to evaluate this term for an r = constant surface with r > 0, which corresponds to a timelike 
surface, and then take the r —^ 0 limit. 

We shall first assume that the variations preserve the GNC form, i.e, we will only vary the functions 
present in the GNC metric, viz. (a,/3^,q"^^). Further since the metric only has the combination ra 
and r^A occurring in it, the variations rSa,r6l3A will vanish in the r —>■ 0 limit. So we are essentially 
restricting ourselves to variations with S£°' = 0, 0 at this stage. (We will describe a more general 

situation later on.) 

When the evaluation of is carried out and the limit r —>■ 0 is taken carefully, we get the result 

- hab6{V\h\iK^^ - Kh-^)) = duiVq^-^) + Vq[-‘^Sa - ^Sg^^dugAB - q^^du{SqAB)], (179) 

where q is the determinant of the 2-metric qAB- The first term in Eq. (179) precisely cancels the 
^/\h\Da{6n‘]_) term so that we readily obtain 

QVWI = v^[-2(5a - ^Sq^^duqAB - q^^du{SqAB)] (180) 

This result may be rewritten in the following form: 

Q^/\h\ = -2y/q5a - duSqAB - y/qdu[Silny/q)], (181) 

We can also derive the same result starting from Eq. (175). Under the variations that we are considering, 
each of the terms in Eq. (175) have the following limits as r —>■ 0: 


V\h\Da (Snl) 
S (27f 
- Kh,,)5h^^ 


-du ( yg— 


5 ( 2duy/q + 2ay/q - y/q 


dua 


^-^duy/^+^y/^5q^°d^qCD 
P^5y/qduy/q + 2.5y/q - 


1 dua 

2 a 


The difference of the last two terms can be simplified and written as 


5{2K^\) - - Kh,j)5h^^ 


5a\ 
a ) 

du6qAB - y/qdu (Slny/q) 


-2Say/q + du y/q 


(182) 

(183) 


(184) 


(185) 


Combining with the total derivative term, we again arrive at Eq. (181) as the null limit for Eq. (175): 


qVW\ 


s/\h\Da{5nl) + 5{2Ky/\h\) " VW\{K.j - 

-2y/qa - ^y/qq^^duSqAB - y/qdu {Sluy/q) 


(186) 
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For the Gaussian null coordinates, we choose the normal Va to be £a = dar^ the surface gradient to the 
r = 0 surface. Then, the normal vector has the following elements in (u,r,x^,x^) coordinates: 

4 = (0,1,0,0); 4=5“" (187) 

(Note that this is not the £ we used in Appendix C.l in constructing the GNC metric.) The auxiliary 
vector (see Appendix A.3.2), defined by the two conditions iak‘^ = £°‘ka = —1 and kak°‘ = 0 on the null 
surface, can be chosen to be ka = —Qar and = (0,—1,0,0). Then, we shall have the induced metric 
on the null surface defined as, qab = gab + £akb + ibka- The explicit form of the induced metric on the 

T* - II IQ 

qab = qABStSl^, q<^^ = q^^5\5l, (188) 

while = (img(0, 0,1,1). Next, let us look at the second fundamental form for the null surface (see 
Appendix A.3.5): 

0a6 = (189) 

Calculating 0ab for our case, we find that the non-zero components are 

&AB = ^duqAB, (190) 

so that the contraction of the above tensor Oab leads to the the trace 

0 = ^q'^^duqAB ■ (191) 

Thus, we get the following expression for (0ab — Qqab)Sq°'^ in GNC: 

{Qab - Qqab)5q“^ = IrduqABSq^^ + -du{^)6{yfq) (192) 

Z q 

Comparing this result with Eq. (181), we find that, for GNC, the surface term is expressible in terms 
of the induced quantities Qab and 0 in the following fashion: 

SAnuii = J ducfx[-2y/qSa-d{2y/qQ) + {Qab-Qqab)Sq‘'^^] 

= J du(fx[-S[2^/q{Q + a)] + [Qab - {Q + a)qAB]Sq'^^^/q] (193) 

Therefore, for the variations we have considered, we only need to fix qAB on the surface, by con¬ 
struction. Even though we varied the six components in the GNC metric, a and appear in the 
combination ra and r/3^. As the null surface remains r = 0 even after variation, when we take the limit 
r —>■ 0 contributions from 5a and 5j3^ vanish. Since the components of 1°“ is given by Eq. (187), this 
corresponds to variations with 5£°^ = 0 and our general result in Eq. (37) shows that we only need to fix 
qAB on the surface. So everything is consistent. In order to get the full structure of Eq. (37), we need 
to consider unconstrained variations of the GNC line element, which we shall take up next. 

C.1.3 Surface Term in GNC for Unconstrained Variations 

We shall now consider the null surface described in the GNC coordinates but allow the variations to be 
arbitrary but finite. That is, the variations are not restricted to be of GNC form. We shall start from 
the expression for the surface term in terms of the variation of the Christoffel symbols given in Eq. (136), 
Appendix B: 

V^gQivc] = V^gva{g^^5Tl, - . (194) 

Substituting the base metric as the GNC metric and Vc = dcV, we obtain 

^gQ = ^ (5r4 + q^^S^AB - - STi^) (195) 
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Evaluating the Fs for arbitrary variations from the GNC metric, we obtain the following expansion: 




driSguu) - -^duqABSq - q duiSqAB) 


Sq^ 


q ^rqAB dui^^Sur) 


+ 5/“ -2a - 


-,AB 


-duqAB 


-Pc - ■^q'^^dcqAB + q^^dAqBC 


Sq-^ 

q^^dA{6guB)] 


(196) 


It will be fruitful to convert the variations of the inverse metric into variations of the metric. For example, 

= -g^^^g^^Sgab := = Sguu • (197) 

Similarly, we can obtain 

<55™ := - Sqru (198) 

Sq''^ ■■= - q^^SguA (199) 

Using these expressions in Eq. (196) and manipulating the terms, we can obtain 

V^Q =dA {^/qq^^SquB) - ^u{.^/q5qur) 

^rQuu )+25{g 

ur^u y/q) “b cx{qru 1)] 

+ 5quudr{y/q) + yfq SquB + cty/q q^^SqAB - ^PuqABSq^^ - ^/q q^^du{5qAB) ■ ( 200 ) 

We will now show that this result matches with the one in Eq. (37). To do this, let us start with the 
following expression: 


0 + « = V^r + 


= V^r + -k<^daq^ 


( 201 ) 


Varying it, we arrive at: 


25 (0 + k) = 2(5 (V„r) + Sk’^daq"''' + k’^daSq''^ 

= 2duSq'^^ + 9^(55™ + 2dA5g^^ + 2Sq^^du In ^ + 2Sq^^dr In 7=^ 

+ 25g^^dA In y/—q + 28^5 In \/—q — 2a5q'^^ 

— 28u5qur 8r5guu 2cx5qur 2/3 5quA 28 a ^guB^) ^gur (5 8uqAB') 

- Sguu {q'^^8rqAB) - q^^SquB {q^^8AqcD) + 2a„(51nv^ (202) 


Using this result, the variation of the counter-term becomes 

2 ^ [\/^ (0 + k)] = 2 y /^5 {Q + k) + 2 -v/^ (0 -I- k) dqur - y/^ (© + k) qABSq^^ 

= yf^ - ‘ 2 ‘ 8 uSqur - drSquu - 28 a {q^^SguB) - Sgur ( 2 a -b 20 -b 2 k -b q^^ 8 uqAB) 
- Sguuq"^^8rqAB - SguB (2/3^ -b q^^q^^8AqcD) + 28u5 In - (0 -b k) qAB^q^^ 


(203) 
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Adding to Eq. (200), we get: 


+ 2S [v^(0 + k)] = -iy/^duSgur - 9a {q^^Sgus) - ^gur {^q^^duqAB + 20 + 2k^ 

- 5guu (^q^^drqA^ — P^SguB — -^q^^SguBq^^dAqco 

+ 2duS In ^J—g — (0 + k) qAB^q^^ — 0ab<^9^^ ~ q'^^du^qAB 
= [0AB - (0 + k) qAB] Sq^^ - 2y/^du5gur - V^Oa (q^^SguB) 

+ SguB + ^q'^^q^^dAqcD^ 

- 2duS In - du {q^^SqAB) 

= [0AB - (0 + k) qAB] Sq"^^ - du (V^^gur) “ 9a {s/^q^^SguB) 

- 2v^ (0 + k) 5gur - V^Sguu (^^q^^drqAB^ - V^P^huB (204) 

In this expression the total divergence term, the counter-term involving 0 -|- k and one involving the 
variation with [0“^ — (0 -I- /c)g“^] are easy to interpret and we have seen them earlier. The extra term 
which appears, in addition to these, has the following form: 

Qextra = +2\/^ (0 + k) (5^“’' drqAB^ + V^^B^g''^ (205) 

These terms have the variation of 6i°' = Sg°'^. One can also show that the coefficient of matches with 
the one in our general expression derived earlier in Eq. (37). Hence the above result can be taken as a 
verification of our general result. If we had restricted our variations to GNC parameters only, then we 
would have Sg^^ = = 0, 5g'^‘^ = 2r5a = 0 and Sg"^^ = 0. Hence we would retrieve our earlier result 

that only needs to be fixed. 


- Sgur (30 -I- 2k) + 5g, 


uu I 2^ drqAB 


C.2 Null Surface Foliation (NSF) 

In this section, we introduce another parametrization for the spacetime metric near a null surface. The 
fiducial null surface is now taken to be a member of a set of null surfaces, unlike the case for GNC where 
only the r = 0 surface in the family of r =constant surfaces was taken to be null. We shall call the 
resulting form of the metric the null surface foliation (NSF) metric. Here we shall give the main outline 
of its construction. The details can be found in [39, 40]. 

C.2.1 Construction 

Consider a four-dimensional spacetime manifold specified by [M,gab)- Let our fiducial null surface be 
one member of a family of null hypersurfaces in this spacetime. Using the fact that null surfaces are 
spanned by null geodesics, it is possible to introduce a natural coordinate system adapted to our family 
of null hypersurfaces in the following manner. We first select one of the co-ordinates, say , such 
that = constant represents the set of null surfaces with a choice of = 0 on S, our fiducial null 
surface. Then, we choose a spacelike Cauchy surface, Ej, and denote the intersection of S with Et, a 
two dimensional surface, by St- Let us define the coordinates on this surface as x^ = (x^,a:^), with the 
two corresponding basis vectors lying on St denoted by ba = djdxA (with the foreknowledge that the 
fourth coordinate will be chosen to be constant on St). At every point P on St, there are exactly two 
future-pointing null directions orthogonal to it, among which one direction will lie on S. Let £ denote a 
null vector field tangent to this direction at every point on St- Thus, i will satisfy the relation I.ba = 0. 
The same exercise can be repeated replacing S with every other 0:3 = constant surface so that we will 
have the coordinates, xi and X 2 (as well as 0 : 3 , which has already been defined all over spacetime), and 
the vector field £ defined all over Ej. Within this arrangement, we can construct a coordinate system 
near E* in the following way: (a) Take the coordinates x^,x'^,x^ as constant all along the null geodesics 
passing through every point P{x^,x^, x^) on E* and moving in the direction of and (b) take t to be an 
afSne parameter distance along these geodesics with, say, t = 0 on Ej. Let us order our four coordinates 
as {t,x^,x'^,x^). In this coordinate system, the null vector is given by £“ = (1,0, 0,0). Thus, the null 
condition = 0 implies gtt = 0- On the other hand, the geodesic condition with affine parametrization. 
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= 0, gives g°‘'^dtgtb = 0; i.e. dtgtb = 0. In other words, gtb remains constant along the null 
geodesics. Since £ ■ ba = 0 on the initial surface, we must have gtA = 0 everywhere. Thus, the form of 
the line element near the null surface may be written as 




Mdx^ 

N 


eNdt 


dAB (dx 


^dx^) (da 


^dx^) , 


(206) 


with e = ±1. The t = constant surfaces are taken as a spacelike surfaces. The a;^ =constant surfaces 
are null with the metric given by = qAsdx^dx^. The inverse metric and Christoffel symbols of this 
metric are given in Appendix D.3. 

This metric may be rewritten in the following condensed form: 

ds^ = M^{dx^)^ + 2Medtdx^ + 2qAB'm,^dx^ dx^ + qABdx^dx^, (207) 

with /N"^ + qABBX'^m^. This metric has seven parameters, one more than the GNC metric in 

Eq. (170). Since it is possible to eliminate four metric coefficients out of the ten independent coefficients 
required to parametrize a symmetric four-dimensional metric by using the four coordinate transforma¬ 
tions, there must be a coordinate transformation that we could make to reduce the number of parameters 
to six. This coordinate transformation has been carried out in Appendix D.2 leading to 

ds^ = (M^ — 2 ef ){dx^)'^ + 2 edldx^ + 2 [qAB'm^ — et )dx^dx^ + qABdx^dx^ 

ox^ ox^ 

= M^{dx^)'^+ 2 edtdx^+ 2 qABm'^dx^dx^ + qABdx^dx^ (208) 

Another way to get rid of M is to go over the construction of the coordinate system once again and 
note that the condition dtga = 0 out of the four conditions dtgta = 0 has not been used. Now, once 
the coordinate system (a:i, X 2 , x^) has been defined on Et and we have decided that Et is a t = constant 
surface for the fourth coordinate t, we have the basis vector 63 = dfdxs corresponding to the X 3 
coordinate at every point on Et. We can then demand i.e^ = 1 everywhere on Et. The four conditions, 
£.£ = 0, £.63 = 1 and £. 6 ^ = 0, fixes £ uniquely. Then, we shall have M = gts = £.63 = 1 on Et and, by 
virtue of dtgts = 0, everywhere else. Thus, from now on, we will assume that the parameter M = 1 and 
do the rest of the calculations. The final form of the metric we shall use is 

/dr^ \2 

ds^ =—N^dt^ + + eNdtj + qAB (dx"^ + m^dx^) (dx^ + dx^) , (209) 

which on expansion becomes 

ds'^ = M'^{dx^Y + “^Medtdx^ + 2qAB'm"^ dx^ dx^ + qABdx^dx^, ( 210 ) 

with = l/N"^ + qAB'm^m^. 


C.2.2 Surface Term in the NSF Metric under constrained variation 

In this case, with the coordinates {t,x^,x^,x^), we shall take the normal to the null surface as 

£a = 9^X3 = (0,0,0,1); r = (e, 0,0,0). (211) 


We shall choose the auxiliary vector 

m2 j^2 

= (-^,0,0,-l); ka = (-e,-mi,-m 2 ,- —), 

where we have dehned M^ = (-^ -I- rn?). One can verify that k°‘ka = 0 and £“fco 


Aj^B 


qabdx°‘dx° = qABdx^dx 


and 

where q"^^ 


q°'^dxadxb = m^dxQ — 2em"^dxodxA + q"^^ dxAdxs 
is the inverse of qAB- Finally, we have the mixed form we have 

qbdxadx^ = —emBdxodx^ -I- SgdxAdx^ 


( 212 ) 

— 1. Then, we have 

(213) 

(214) 

(215) 
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The non-zero components of 0ab are 


(216) 


and 



(217) 


where q is the determinant of the 2-metric qab ■ We shall now write down the surface term. The surface 
term, y/—g can be decomposed using Eq. (16). The different terms in this expansion give: 


= -dt (218) 

-<5 (2^^ V.r) = -(5 (2MVg 0) = [2dt (e^)] (219) 

( 220 ) 

-V^gabVcrSg^^^ = -v ^0 [2eM5 (^) + qabSq^^] ( 221 ) 


Our aim was to see if we can fix just on the surface, just as we only need to fix h°‘^ on the surface 
for a spacelike or timelike surface. We have variations of M in Eq. (221). But we have already seen in 
Appendix C.2.1 that M can be set to unity with a suitable choice of coordinates. Then, the surface term 
takes the form 


SAnuU= [ dX(fz±{-S[2y/q{Q + K)] + y/q[0AB-i0 + K)qAB]Sq'^^} ■ (222) 

JdV 

We see that we get back the exact expression we had obtained for GNC in Eq. (193), with a and qAB 
replaced by their corresponding quantities k and qAB- Note that the variations we have considered 
sets 5(°' = 0 and hence we pick up only the remaining terms obtained earlier. To reproduce the terms 
involving 5^‘^ we need to consider a more general form of the variation, which we shall now take up. 


C.2.3 Surface Term in NSF for unconstrained variation 

We shall now take the original metric to be the NSF metric but allow the variations to be arbitrary, i.e. 
the variations are not restricted to those which preserve the original NSF form. We shall start from the 
expression for the surface term in terms of the variation of the Christoffel symbols given in Eq. (136) in 
Appendix B: 

^gQ[v,] = ,/^gv,{g^^5Vl, - . (223) 

Substituting the base metric as the GNC metric and Vc = dcX^, we obtain: 

QH = g^*5Tl + + 2g*HTl, - g^^5T\, - g^Hvi, - g^^dVl + q^^ (224) 

Let us now evaluate each term individually and then put them together. For that, we have the following 
individual expressions: 






leg^^Sgtt 


-em^ {Sg^^dtruA + Sg^^dtqAB + edASgtt) 


Sg-^-^dt ( m + j;p] +^9 dtniA + edsSgu 


,AB 


- edtSqAB + (^5 {-daqAB + dArriB + dBrriA) 


Sg'^^ {-dcqAB + dAqcB + dBqcA) + <^5 {-dtqAB) + e (dASgtB + dBSgtA) 


-ppdtSgtt - edsSgtt - em^ {-dASgtt + 2dt5gAt) 


i [Sg^’^dtruA + Sq^^dtqAB - em^dASgtt + q'^^dtSqAB] 
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Then, substitution of these results in (Eq. (224)) leads to the following expression; 
Qbc] = ^g^^dtSgtt - ern^ {Sg^^dtmA + Sg^^dtqAB + edASgtt) 


e 

2 

1 

2 ^^ 


Sg^^dt + ^] + Sg^^dtruA + ed^Sgu 


AB 


- edtSqAB + 5g'^'^ {-daQAB + dAms + dsmA) + 5g^^ {-dcqAB + dAqcB + dBqcA) 


+ ^9^* {~9 ^dtQAs) + € {dA^gtB + dnSgtA) 


-jpdtSgtt - edsSgtt - {-dA^gu + “IdiSgAt) 


- 2 [^ydtmA + Sq^^dtqAB - em^dA^gu + q^^dtSqAB] 

Now the variation of the expansion scalar can be obtained as 

2(50 = 2(5 {Vat) 

= 2dtSg*^ + 2dz5g^^ + 2dASg^^ + 25g°‘^da In ^/^ + 2dt {S In ^/^) , 
which can be used subsequently to obtain the following result: 

2 ^ (V^ 0 ) = ‘2\A^dt5g*^ + 2y^d36g^^ + 2y/^dASg^^ + 2dg°-^da\t^ 
+ 2A^dt {6\n^/^) - 2^/^QSg*^ - A^^qAB^q"^^ ■ 

Also, we have the results 

<5^“ = SgtA = -egtASg^* - eqABSg^^ - egAs^g^^ ■ 

Thus, we arrive at the following result: 

V^Q + ‘2S = |9t - dA {yf^m^Sgtt) 

+ edA {y/^q^^SgtB) + (©as - ©QAb) Sq^^ + yf^Qextra 


(225) 

(226) 

(227) 

(228) 

(229) 


As usual, we have the standard counter-term and variation proportional to Sq^^ which are by now 
familiar. In addition to boundary terms and counter-terms, we have the following extra terms; 

V^Qextra = -y/^Sg^^q'^^ dtqAB - €y/^Sg^^ [mAq^^dtqcD + dtniA] 


yA^5g 
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AB 


dtqAB 


(230) 


These terms are all proportional to i.e. to St, as to be expected, with the proper coefficients dictated 
by Eq. (37). This matches with our earlier conclusion. If we restrict to variations of NSF parameters, 
then it is clear that St = 0. Then, in accordance with the general result presented in Eq. (37), only 
variations of need to be fixed on the boundary, as we have shown explicitly in the previous section. 

D More on GNC and NSF 

D.l Inverse Metric and Christoffel Symbols in GNC 

With coordinates ordered as {u,r,xA, the metric as presented in Eq. (170) and its inverse may be 
written in matrix notation as 


—2ra 1 —tPa 
gab=\ I 0 0 

-rjSA 0 qAB 


0 1 0 
g°-^ = I 1 2ra + 


0 r/3"' 


^AB 


(231) 


where q^^ = (g“^)"^^ is the inverse matrix of qAB, /3^ = q^^P b and /3^ = Pa(S^- The determinant of 
the metric is g = q, where q is the determinant of the 2-metric qAB- Let Da be the covariant derivative 
operator associated with the two-dimensional metric qAB ■ For example. 


DaISb = SaPb — 


(232) 


37 










where 


{—daqAB + dAQBD + duquA) 


^ AB — 


^CD 


2 


(233) 


and DaP = PaP for the two-dimensional scalar p. Introducing a = —2ra and Pa = —tPa, the Christoffel 
symbols have the following expressions [38]: 


Ku = -\dra 

Ka = -\drPA 

Tab = —^^rqAB 

-pli Y'U Y^u n 

^ ur ^ rr ^ rA ^ 

r;;, = lidra-P^drPc) 

= lidrh-P^drqCA) 

^uu = -]^{P^Pc-a)dra+]^dua+^P^dca-P^duPc 

^'ab ~ ~2 {P^Pc ~ d) drqAB^ + 2 {J^aPb + DbPa^ 

^uA = -^{P^Pc-a)drPA + ^dAa-^P^ (duqAB+dAPB-dBpA) 
= 0 

r^c = Ip^drqBc + ric 

^Bu = ^P^drPB + ^q^^duqBC + ^q^^ {dBPc “ PcPb) 

= ^q^^drqBC 

^uu = ^P^dra - ^q^^dca + q^^duPc 

^ur = \q^^drPc 

= 0 


(234) 

(235) 

(236) 

(237) 

(238) 

(239) 

(240) 

(241) 

(242) 

(243) 

(244) 

(245) 

(246) 

(247) 

(248) 

(249) 


D.2 Eliminating the Extra Parameter M in the NSF Metric 

The NSF metric was given in Eq. (207) and is of the form 

ds^ = M'^{dx^Y + ‘^Medtdx^ + 2qAB'm^dx^ dx^ + qABdx^dx^, (250) 

The condition that the null vector i°‘ = d/dt is affinely parametrized leads to dtM = 0. So, we have 
M = M{x^,x‘^, x^). 

Let us make the coordinate transformation to a new coordinate i = tM{x^, x^, x^) so that 

DM 1 DM y dM n , , 

Mdt = dt — t-^—dx^ — t-^^dx^ — t—^dx^ . (251) 

ox^ ox^ ox-^ 

The metric in the new coordinates would take the form 

ds^ = {M^ — 2et ^ g ){dx^Y + ‘2^f-dtdx^ + 2{qAB'm,'^ — et ——^)dx^dx^ + qABdx^dx^ . (252) 

We can see that this metric is of the same form as we would have obtained by putting M = 1 in Eq. (250). 
P/Pt is a null geodesic with affine parametrization, as we would have expected since t is a linear function 
of t and hence also an affine parameter. Thus, we have managed to encode the same information as 
in Eq. (250) with one less parameter. But one difference with Eq. (250) is that while P/Px^ was taken 
to be a spacelike vector in Eq. (250), the vector P/Pxs in Eq. (208) (which is a different vector) is not 
necessarily spacelike. 
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D.3 Inverse Metric and Christoffel Symbols for the NSF Metric 


The NSF metric corresponding to the line element in Eq. (207) and its inverse in matrix form are given 
below in the coordinates ,x^): 


/ 0 0 eM 

gab= [ 0 qAB rUA 

\ eM ruA ^ 



(253) 


where = {q~^)^^ is the inverse matrix of qAB^ = q^^niB and = mAm^. The determinant 
of the metric \s g = M^q, where q is the determinant of the 2-metric qAB- 

Introducing Da as the covariant derivative operator corresponding to the two-dimensional metric 
qAB-, as in Appendix D.l, the Christoffel symbols for the above metric are 


r* 

^ At 

r* 

^ AB 

pt 
^ X3t 


pi 

^ X 3 A 


Bt 


BC 


tX3 


Bx3 


pa;3 

2:3X3 


T-'2;3 

^ Ax3 
^ AB 


p2:3 
^ tX3 


dtlnM 

i^AlnM -t -^qABdtm^ 


^dtqAB + ^ + DBmA) - ^d^qAB 


2Af2 

1 


-m^dA In M -h 1 irPT I + TTlTJ^Adtru 

M2 


M2 


2M 


i9a - 




iV2 

e 


2M 

{OaM - edtiriA) 


[dATTlB - OBmA) - 


2M 

1 


2M 


m2 


2M 


dAirP + 


M2 

Ip 




ATS 


2M 


m^m^dsqAB 


2^ 


AC 


dtqBC 


em 


dtqBC + ^BC 


—2^ 


AB 


-,AB 


ObM ' 


-^asM 

0 


em 


-dtmB 


2M 


dtmB 


M2\ 


em 

- Of I m~ -h „ , 

2M \ m) 

- a 1 ^ ^ 

“TTlTT^t -I- 


2^ 


AC 


{dBmc - dcmB) 


1 


^ 1 - In M - ^q^^dB q^'^dsmB 


2^ 


m2 


AC 


d^qBC 


iV2 




2M 


7V2 


ed^ In M 


-^dtqAB 

= r^l = 0 . 


(254) 

(255) 

(256) 

(257) 


(258) 

(259) 

(260) 
(261) 
(262) 

(263) 

(264) 

(265) 

(266) 

(267) 

(268) 
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E Interpretation of the Counter-term in the Non-null Case 

We have seen that the counter-term in the null case can be interpreted as the difference in the 2-surface 
areas orthogonal to £ and k at the boundaries of the null congruence on the null surface (see Eq. (45)) 
when = 0 everywhere. The counter-term to be added to the boundary in the non-null case is the 
integral over the boundary of the expression 2A/|7i|Von“ (see Eq. (164)). We shall work with the case 
where the boundary is a spacelike surface, a t-constant surface in coordinates (t, j/i, 2 / 2 , 2 / 3 ), and provide 
an interpretation for the counter-term. Working along the same lines, a similar interpretation can be 
given for the counter-term on a general non-null surface. The metric has the ADM form [21, 1] 

ds^ = -N^dt^+ h^is{dx°+N°dt){dx<^+N^dt) , (270) 
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with the inverse metric components given by 

1 


5 “ = 


Ar2’ 


9 iV2 ’ ^ iV2 


(271) 


The normal according to our conventions (see Appendix B) is Ua = dat. The counter-term integrand 
VhVa n°‘ can be manipulated as follows: 


VhVan‘^ = 


y/h 

7^9 


da = -^da [Ny/htf-^ = n°‘da + ^da {Nn°-) . (272) 


We have n“ = = (1/A^, N^^/N), the second term in Eq. (272) can be written in the following form: 


—da (Nn^) = —daN° 
N ^ ’n 


(273) 


Let e“, a = 1,2,3, represent the coordinate basis vectors on the 3-surfaces normal to n“. They satisfy 
e“na = 0. We shall consider the case where these basis vectors are Lie transported along So we have 

n^dbe't - = 0 . 


K,e“]=0 


(274) 


In the coordinates where e“ are basis vectors, we would have e“ = (5“. Hence, the above condition 
reduces to 


daTi^ = 0 => daN = 0 and daN°‘ = 0 . 


(275) 


These conditions also imply that n“ are tangent vectors to affinely parametrized geodesics. Using Ua = 
NS/at and S/oN^bt = VbS/at, we can reduce the geodesic equation to the condition 

h'^^dahiN = 0 . 


n°N/anb = 0 


(276) 


Since has only the spatial components, Eq. (275) implies the RHS is zero. Hence, n“ should 
satisfy the geodesic condition. Using Eq. (273) and Eq. (275), Eq. (272) reduces to 


'/hVan°‘ = (^) ’ 


(277) 


which is the change in the volume of a 3-surface element along Taking r to be the parameter along 
the integral curves of n“, we may also write the above equation as 


y/hS/an°' = 


5(v^) 

dr 


(278) 


Unlike the null case in Eq. (45), where the derivative was on the surface, this derivative is off the boundary 
surface and hence will not get integrated in the boundary integration to be interpreted as the difference 
in the volume of the 3-surface element at two different points. 


F Details of Various Calculations 

F.l Derivation of Eq. (49) 

The boundary term in the non-null case has the following expression: 

y/\h\Da {SnD = da 

= da f y/\h\ + 2rSN) 

= da {Sr + 2rd In N)] 

where we have used Sia = 0 to get to the second line and Eq. (145) in the last line. On the other hand, 
the surface term on the null surface in Eq. (37) can be manipulated as follows: 

da = da [V^n^bK] = da + da 

= da + da [V^k^^ibSi’’] 

= da 
ab i 


(279) 


(280) 


where we have used 5ia = 0 to write Si°[ = 5£^ -I- g°^°5£a = 5£°' and used 5(P) = £a^£^ = 0 on the null 
surface. 
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F.2 Derivation of Eq. (50) 

The variation of the counter-term can be written as: 

<5 = -,5 (2v17^Van“) 


= -25 


= -25 


{^99^*N) 


= -25 [da - 2<5 [da (IniV) 

= -2,5 [da (^r)] - 25 [^tda (IniV)] 


= -25 [da (y^r)] + 5 




(281) 


where the last line has been obtained using the relation g'^'^ = e/N^. In the limit iV —>■ oo, it is clear 
that all terms in Eq. (281) are finite under our assumptions. For a null surface, we have the following 
result using Eq. (128): 


5a (y=^r) = V^v.r = I (0 + k) - I (0 + «) _ 


(282) 


As we are interested in considering variations such that is kept zero on the surface (since we demanded 
5{(.ai°‘) = 0 in obtaining Eq. (37)), the last term in Eq. (282) can be manipulated as follows in the null 
limit: 


2 2 


= 


(283) 


where we have used k°‘ia = k'^ = —1 in the last line. Thus, we can write the null limit of Eq. (281) as 




5 {2^\Kj -25 (0 + k)] - 5 [V^d^g^^] + 5 




— —26 [v ^—g (0 + + (5 




(284) 


= -25 (0 -h k)] - 25 [y/^ e°‘dc In N] 

F.3 Derivation of Eq. (51) 

We will first calculate the quantity hij5h^^: 

hij5h^^ = [hij5g^^ — 2ehijn^5n^) 

= 

= (gij^g^^ - en,nj5g"^) 

= {g,,5g^^ - tNHdj5g^^) 

= {g,j5g^^ - tNHg**) 

= {g,j5g^^+25\nN) 

On the null side, we have the following result: 

gab5g“*' = qab5g°‘'' = gab5g°‘'" + 2lakb5g°‘^ 

= gabSg'^^ + 2kb5g’>^ = gab5g‘^^ + 2kb5t 

Now, if we look at the derivation of Eq. (50), it is easy to see that the expression is valid without the 5 
too. Thus, we have 

V^iC’'^°-^[(0 + K)+r5clnfV] . (287) 

We are ready now to attack the last term in Eq. (46): 

^/\h\Kh,j5h^^ -v/^[(0 + k)+ In TV] x [g^j5g^^ + 251n Af] 

= (0 + ^) qab5q°'^ + 2-/^ (0 -I- k) kb5& — 2y/^ (0 -|- k) 5 In TV 

-V^^“5„lnTVgy55*^ -28/^(r5clnTV) 51nAf (288) 


(285) 


(286) 
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F.4 Derivation of Eq. (52) 

We have: 


- ^/\h\K,jSh^^ = -^/\h\KiJSg^^ +2e^/\h\Kijn^Sn^ 

= -^/mSg"-^ [-ViUj + eriin^Vrnnj] 

= [^^^3 + In iV - - eN^W'^V„^ In N] 6g^^ . (289) 

The third term in Eq. (289) is 

- eN^eirVm^jSg"^ = 


where we have used the symmetry of Va(-b in the second step. This term adds with the second term in 
Eq. (289), while the last term in Eq. (289) gives 

-eN^tdjr'{y^\xiN)5g^^ = -eN^l^r {V ^\uN) 5^^ 

= 2r{d^\nN)5\nN , (291) 

where we have used (.a5t°' = 51^^ = Sg'^'^ in the first step. Adding everything up, we obtain 

-^/\h\Kij5h^^ = ^ [V^^j5g^^ + 2 (9, In N) St + 2^ (9, In N) (5 In7V] . (292) 

All terms here are finite in the null limit. In the last term, SlnN is finite under our assumptions and 
£'‘ {di In A^) has been shown to be finite in Eq. (281). {dt In TV) 5t can also be easily shown to be finite. 
The first term in Eq. (292) needs to be decomposed in the null limit. This is done as follows: 


{djg’*”*’) 5& 


2g<l><l> 

SfdjlnN 


(290) 


\/iijSg^^ '■=° Vdj [Sq"^ - 2S (tk^)] 

= - Sk^dd"^ - 25ek^Vd, 

= OijSq^^ + 25(.d^Vikj . 


(293) 


In the last step, we have used Eq. (33) for the first term and 5£a = 0 was used to put Sk^dd"^ = 0. 
Putting it all together, we arrive at 


- y/\h\K^jSh^^ = Q^JSq^^ + 2StdVdj + 2 (9^ In TV) 5e + 2f {d^ In TV) 5 In N 


(294) 


F.5 Derivation of Eq. (54) 

The extra terms in Eq. (53) are 


Extra Terms = 2da [-v/— gP^iiln TV] — 25 \\/~9 ^°‘9a InTV] — 2y/—g (0 + d) (5 In TV 
- yf^ (^“3c In TV) - 2yf^ (^“5^ In TV) <5 In TV 

+ y/^h {di In TV) 5t + 2t {di In TV) 5 In A^ 


(295) 


We shall start by working on the first, second and the fourth of these extra terms. These are 

2d^ [v^r 5 In N] - 25 [y/^ In A^] - y/^ (£“5^ In N) g,j5g^^ 

= 2[da{y/^r)5lnN - y/^5rdalnN] (296) 

Substituting back, we have the extra terms as 

Extra Terms = 2{[d^ (\/=5^“) - \/=5 (© + «)] 5 In A^ + y/^5i‘l>d^ In A^ + yf^ {£^d^ In A^) ,5 In TV} . 

(297) 
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Now, from Eq. (282) and Eq. (283), we know that the first term in Eq. (297) can be manipulated as 


[do. - ^(0 + «)] 5\nN= -e-n)-d^ 'Jin TV 

= 6lnN 

= Jin TV 

= (y=5^^,^lniv) JlnTV 

= In N) 6 In N, (298) 

where we have made use of Eq. (142) and £“ = g“'J. In the second step, we have put the term 
— {^/—gi^d^y/—g) (5 In TV to zero as we know that = i°‘^a = 0 on the null surface while \/—g and 
its derivatives are assumed to be finite everywhere and d In TV is finite under our assumption about the 
behaviour of g'^^ near the surface. 

Thus, the first and the last terms in Eq. (297) are identical. The middle term should also be then of 
the same form with an extra factor of —2. Indeed, we have 


In TV = — {d^ In TV) ,5 In TV = {d^ In TV) d In TV , 


so that Eq. (297) reduces to 


Extra Terms = 0! 


(299) 


G Working with a General ia = AVa(t> 

In this appendix, we shall consider the general case of our null normal being of the form ia = = Ava, 

for an arbitrary scalar A which may depend on the metric. We no longer have the results Va^b = Vb^a 
everywhere and 6ia = 0 everywhere that we had used profusely in Section 3.3. But we do have the result 

Sia = SAVafj) = S In A ia (300) 

From Eq. (18), the boundary term on the null surface is 

= ^Q[L] = J {v^Vc[J^l] - 2S{V^Van + yTgiVait - gat^dlSg'^A ■ (301) 

Labelling the first term as ^—gRi , we have 

V^Rl = ^^a[6il] = jdaiV^Sil] = dal^Sil] - d^Si^da (302) 

We shall now use the projector 11“^, to separate out the surface derivatives in the first term. 
da[^6il] = da[^Il%Si\^] - dai^kdbSA] 

= da[^Il%Si\_] -SiLndb [6{ian] 

= da[^Il%6i\^] - ^k^db [J iiat)] , (303) 

where the last step was obtained by using our assumption 6 {iad) = 0 on the null surface. Using this 
expression, we have 

yTgRi = da [^n“bJ^i] - ^k»db [J {iad)] - V^Si^da (304) 

The first term in Eq. (304) is a surface derivative on the null surface as = 0- Th® second term in 

Eq. (304) has variations of the derivatives of the metric. We shall take out the 5 to obtain 

- k^db [J {iat)] = -5 [k^db {iat)] + Sk'^db {iat) . (305) 
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Substituting in Eq. (304), we obtain 


V^Ri =da[^n\dii] - ^6 [k^db (4r)] + ^Sk^db {Ln - V^seida 


=da[^Tl\6£’i] - S 


:^k^db{ian 


V^9 

2A 


[k^db {Ln] + V^k^db {Ln s 




A 


6k^db {Ln - V^seida 


1 
A 

(306) 


Here, all the variations of the derivatives of the metric are in the first two terms, assuming A does not 
depend on the derivatives of the metric. The second term in Eq. (301) is 

v/^i?2 = = -2<5(^VX) + 2y^VarS (307) 

Substituting Eq. (306) and Eq. (307) back in Eq. (301), and using the relation 

V^r + y5a (hi^) = + k^hVai^ = Il\W, (308) 

in three places, the boundary term on the null surface reduces to 

y/^QlL] =5a[^n“,«i] - 2S{^U\VJ^) + ^(Va4 - gabU^d^ci‘^)Sg‘^^ 

+ 2V^n%Vai^s + ^Sk^db {iaH - ^TgSilda Q) , (309) 


When A = 1, the last three terms vanish and this result reduces to the result in Eq. (29). To see why 
the second-to-last term did not appear for A = 1, note that this term will only have normal derivatives 
as £a-^“ is fixed to zero everywhere on the null surface. Hence, only Sk‘^ contributes in our (^(f),y^ ,y^,y^) 
coordinate system. When A = 1, Sk‘^ = 6 {(.ak°‘) = 0. In the general case, fc“£a = —1 means k'^A = —1, 
which gives 

k* = -llA- 5k* =^-^. (310) 

In Eq. (309), we have succeeded in separating out a total derivative on the surface and a total variation 
to remove all derivatives of the metric. The counter-term to be added in this case is the integral over 
the null surface of 


= 2^ {g\Va(^ - rkbVaH^) = 2^ (0 + k) . (311) 

We shall now do the analysis of what is to be fixed on the null boundary in this case. Since Eq. (17) are 
taken to be valid even on variation, the relations = 0 and q^'^ka = 0 are respected by the variations 
and terms of the form iahSq°'^, £akbSq°'^ and kakbSq^-^ would reduce to zero, just as in Section 3.3. Thus, 
we can simplify gabSg‘^^ as follows: 

gabSg'^’^ = gab [Sq'^^ - S [tk^] - 5 {£^k<^)] = qabSq'^^ + 2 {iakb + £bka) s (rA:*-) = + 2£bka5 (rfc'') 

= qabSq^'^ - 2fca(5r - 2£aSk‘^ = qabSq'^^ - 2ka5£‘^ - 26 In A, (312) 

where we have used Eq. (310) in the last step. 

Next, we shall simplify [Va£b) Sg°'^. We have 

(Va4) Sg'^^ = {Va£b) Sq^’^ - S (rfc^) Va4 - <5 (£^fc“) Va4 

= (Va4) Sq^’^ - 6£‘^k^Va£b - Sk^£‘^Va£b " ^£^k<^Va£b " SkH^Va£b 

XUa 

= {Va£b) Sq^^ - 8tk^ {Va£b + ^b£a) " b^Sk^£b " — 

KUO- 

= {'^a£b) 5q"^ - S£‘^k^ (Va4 + Vb4) -nSlnA- —da£\ (313) 
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where we have used i°‘'S/a^-b = (see Eq. (13)) and iaSk°- = i51n^ (see Eq. (310)). The expression 
(Vo^b) Sq°'^ can be simplified as follows: 

(Va4) = C^b (V™4) Sq-^ = (g™ - - k^la) (%" - ^"^b - fc”4) {VM ^9“" 

= (C - (g," - rkb) {VM Sq^^ 

= (C9b - q^e^ktVmin - qbi^kaVmin) 5q°-^ 

= (oab - _ nq^kj^ Sq^^ 

= QabSq'^^ (314) 


where we used iaSq°'^ = S [q°‘^ia) — q‘^’’Sia = — 9“^4<51nA = 0 to get to the second line, kakb5q°'^ = 0 
to get to the third line and the definition of ©ab and k, (see Appendix A.3.5 and Eq. (13)) to get to the 
fourth line. The final result is obtained using ^aq°‘^ = 0 and the fact that q^dm^'^, a derivative on the 
null surface, is zero since is zero all over the null surface. Next, we shall simplify the last two terms 
in Eq. (313) as follows: 

-KSlnA-^daf = -KS\nA-^d^f =-KS\nA-^^d^e^ 

Ld A Zdi\. 

= -(51nA =-,51nA =-51nA(ft: + K) . (315) 

Here, the first step made use of the fact that I? all over the null surface and hence only the derivative 
along (j) contributes and the second step as well as the second-to-last step used Eq. (310). The last line 
used the definition 2k = —(see Appendix A.3.6). The remaining terms in Eq. (309) are 

^ =-2^^^(0 + K)(51nA (316) 

^^6k^db {tat) = 2^^K(51n A (317) 

J\. J\ 

-V^gSt^da Q) = ^ {5tda\nA + g<^^da\nA5tb) 

= {St da In A + tda In A5 In A) 

= ^^{8tda\nA + {K-k)5\nA) (318) 

Eq. (317) can be obtained using the same manipulations that we performed in Eq. (315). Combining all 
the above results, the boundary term in Eq. (309) becomes 

=da[S^A\5t\] - 25 

+ ^ {2ka (0 + k) - k^ {Vat + Vbla) + da In A) dt . (319) 

The changes from Eq. (36) are as follows: i/—g has been replaced with yJ—gjA everywhere. The da In A 
is extra and 2tVckb5t has been replaced with —at + Vbta) 5t. 


A 


(0 + k) 


A 


[(©ab - (0 + k) qab) 5q' 


.abl 


2^IV\Vat5 
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